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EXECUTIVE  SUMMARY 


The  objective  of  the  research  described  in  this  report  is  the  optimisation  of  wavelength 
conversion  of  laser  radiation.  Wavelength  conversion  of  a  specific  laser  is  a  method  to 
enhance  its  applicability,  since  the  wavelength  of  a  laser  determines  its  utilisation.  In 
particular,  the  wavelength  conversion  of  the  Nd:YAG  laser  from  1.06  pm  to  1-54  pm  is 
attractive  since  the  1.54  pm  is  eye-safe  in  contrast  with  the  eye-hasardous  1.06  pm.  For  the 
armed  forces,  applications  are  found  in  laser  range  finders,  trackers,  and  target  indicators, 
which  can  be  operated  safely  using  1.54  pm  radiation.  Since  atmospheric  transmission 
at  this  wavelength  is  slightly  better  than  for  1.06  pm  and  detection  systems  have  a  higher 
sensitivity  at  this  wavelength,  laser  range  finders  have  an  enhanced  performance  at  1.54  pm. 
Contact  with  Dutch  and  International  companies  has  been  established  for  the  production 
of  a  1.54  pm  laser  range  finder  based  on  the  SRS  technology. 

We  studied  the  technology  of  wavelength  conversion  by  the  stimulated  Raman  scattering 
(SRS)  process.  Using  high  pressure  methane  as  a  SRS  medium,  an  energy  efficiency  of 
45%  was  obtained.  Due  to  the  difference  in  photon  energy,  the  theoretical  maximum  of  the 
energy  efficiency  is  69%.  A  numerical  model  has  been  developed  to  assist  in  the  optimization 
and  to  analyze  the  experimental  results.  Experimental  validation  of  the  numerical  results 
has  been  carried  out  with  success.  With  the  use  of  the  numerical  model,  it  is  expected 
that  other  wavelength  conversions  can  be  optimized  faster  and  less  expensive  than  by  the 
experimental  trial-and-error  approach. 


Conversion  to  the  mid-infrared  wavelength  of  3.6  pm  is  under  w»  using  a  modified  Nd:YAG 
laser  at  1.44  pm  and  converting  this  wavelength  by  SRS  in  hydrogen  gas.  Application  of  a 
3.6  pm  laser  is  found  in  ranging  and  in  jamming  of  sensor  systems  that  operate  in  the  3- 
5  pm  atmospheric  window.  With  the  developed  numerical  model,  the  wavelength  conversion 
to  3.6  pm  will  be  optimised. 
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ABSTRACT  (ONGERUBRICEERD) 

A  three-dimensional  finite  difference  method  is  presented  for  the  propagation  through  a  gain  medium  of 
beams  focused  by  an  astigmatic  lens.  The  method  is  used  for  the  modeling  of  a  Raman  amplifier  that  uses 
astigmatic  lenses.  With  this  method,  die  effect  of  gain  focusing  on  the  amplification  is  studied.  The 
numerical  results  show  a  good  correspondence  with  analytical  models  in  the  range  where  these  models  are 
applicable. 

The  treshold  power  of  stimulated  Raman  scattering  (SRS)  has  been  calculated  for  various  nondifiraction- 
limited  pump  beams  using  the  numerical  model.  It  is  shown  that  the  threshold  power,  for  a  tightly  focused 
beam,  depends  on  die  M2  factor  that  is  related  to  the  pump  beam  quality.  Additionally,  die  influence  of  the 
exact  pump  beam  profile  on  the  threshold  is  small.  Therefore,  die  M2  factor  is  a  useful  parameter  for 
analyzing  the  experimental  threshold.  We  show  that  the  measured  SRS  thresholds  in  methane  correspond 
with  the  calculated  thresholds. 

An  astigmatic  focus  has  been  used  to  increase  the  conversion  efficiency  of  a  Raman  cell  at  high  pump 
energy.  Experimental  and  numerical  results  show  that  die  increased  conversion  is  due  to  die  reduction  of 
cascade  second  order  Stokes.  It  is  shown  that  other  effects,  namely  Brillouin  scattering,  anti-Stokes 
generation,  and  ground-state  depletion  are  negligible  for  our  experimental  setup. 
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SAMENVATTING  (ONGERUBRICEERD) 

Een  drie-dimensionale  numerieke  methode  op  basis  van  finite  differences  voor  de  propagatie  door  een 
versterkend  medium  van  bundels  gefocusseerd  met  een  astigmadsche  lens  wordt  besproken.  De  methode 
wordt  gebruikt  voor  de  modellering  van  een  Raman  cel  met  astigmadsche  lenzen.  Het  effect  van  gain 
focusing  op  de  versterking  wordt  met  het  model  onderzocht  De  numerieke  resultaten  komen  overecn  met 
analytische  modellen  voor  het  gebied  waar  deze  geldig  zijn. 

Het  threshold  vermogen  van  gestimuleerde  Raman  verstrooiing  (SRS)  is  berekend  voor  diverse  niet- 
diffractie-begrensde  pompbuadels  m.b.v.  het  numerieke  model.  Er  wordt  aangetoond  dat  het  threshold 
vermogen  voor  een  gefocusseerde  pompbundel  afhangt  van  de  M2  factor  die  een  maat  is  voor  de 
bundelkwaliteit.  Daarbij  blijkt  de  invloed  van  het  predeze  pompbundel  profiel  op  het  threshold  vermogen 
gering  te  zijn.  Hierdoor  is  de  M2  factor  een  bruikbare  parameter  voor  de  analyse  van  de  experimentele 
threshold.  We  tonen  aan  dat  de  gemeten  SRS  thresholds  in  methaan  overeenkomen  met  de  berekende 
thresholds. 

Een  aatigmatisch  focus  is  gebruikt  om  het  convenie  reodement  van  een  Raman  cel  te  veritogen  bij  hoge 
pompenergie&i.  Experimentele  en  numerieke  resultaten  laten  zien  dat  bet  toegenomen  reodement  te  dankcn 
is  aan  de  afname  van  cascade  tweede  orde  Stokes.  Er  wordt  aangetoond  dat  andere  effecten  als  Brilkmin 
verstrooiing,  anti-Stokes  productie,  en  depletie  van  de  grondtoestand  te  verwaariozen  zijn  voor  ooze 
experimentele  opstelling. 
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1  INTRODUCTION 


Wavelength  conversion  by  stimulated  Raman  scattering  (SRS)  has  been  used  to  generate 
new  wavelengths  from  existing  efficient  lasers  fl].  In  particular,  the  conversion  of  Nd'.YAG 
radiation  at  1.064  fim  to  1.543  pm  by  SRS  in  methane  has  been  studied  at  our  labora¬ 
tory.  The  1.543  pm  radiation  has  the  advantage  over  1.064  pm  that  it  is  eye-aafe  since  the 
radiation  is  strongly  absorbed  by  the  eye’s  interior  before  reaching  the  retina  [2].  Subse¬ 
quently,  we  shall  use  the  terms  pump  beam  and  Stokes  beam  for  the  original  laser  beam 
and  the  wavelength  shifted  beam.  In  this  report,  the  theory  of  SRS  and  the  modeling  of  the 
wavelength  conversion  is  described.  More  information  on  the  experimental  setup  and  ex¬ 
perimental  results  can  be  found  in  the  reports:  FEL-90-B098  (in  Dutch)  and  FEL-92-A125 
13,4). 

Chapter  2  describes  the  numerical  model  that  has  been  developed  at  our  laboratory  to 
study  and  to  optimize  wavelength  conversion  of  laser  radiation  by  SRS.  The  extension  to 
an  astigmatic  focus  is  included  in  the  numerical  model  for  the  analysis  of  our  experiments 
with  astigmatic  lenses.  These  experiments  showed  that  the  use  of  an  astigmatic  lens  or 
a  pair  of  cylindrical  lenses  can  increase  the  conversion  efficiency  of  a  Raman  amplifier  [3]. 
Because  of  the  astigmatic  focus,  the  gain  medium  is  no  longer  rotationally  symmetric,  so  a 
single  radial  coordinate  cannot  be  used.  Therefore,  a  two  dimensional  grid  is  used  for  the 
finite  difference  (FD)  method,  and  this  grid  is  propagated  along  the  optical  axis  resulting 
in  three  space  dimensions.  Previous  studies  have  shown  that  there  is  a  narrowing  of  the 
Stokes  beam  because  of  the  higher  amplification  at  the  center  of  the  pump  beam  [5,6,7]. 
This  effect  is  called  gain  focusing  and  the  analytical  model  describing  it  [6]  is  compared 
with  our  numerical  results. 

Threshold  measurements  and  simulations  are  described  in  chapter  4.  The  SRS  process 
has  no  real  threshold  like  a  laser,  but  an  experimental  threshold  can  be  defined.  Here  the 
threshold  stands  for  the  pump  power  corresponding  to  1  percent  conversion  efficiency.  In 
this  chapter  we  study  the  effect  of  pump  beam  quality  on  the  threshold  power.  By  changing 
the  diaphragm  inside  the  laser,  we  could  change  the  quality  of  the  output  beam.  Because 
the  pump  beam  is  focused  into  the  Raman  medium,  the  SRS  process  will  depend  on  its 
beam  quality  since  the  process  is  driven  by  the  intensity  at  the  focus,  which  is  high  for  a 
good  beam  quality  and  decreases  with  decreasing  beam  quality. 

The  beam  quality  is  related  to  the  Ms  factor,  which  is  a  dimensionless  factor  that  compares 
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the  divergence  of  a  given  beam  with  the  divergence  of  a  (diffraction-limited)  Gaussian  beam. 
The  divergence  of  an  arbitrary  beam  can  be  defined  unambiguously  using  the  variances  a\ 
and  a*  of  the  beam  profile  in  the  two  transverse  directions.  By  comparing  the  rate  of 
increase  of  the  variance  of  an  arbitrary  beam  (i.e.  the  divergence)  with  the  corresponding 
rate  of  a  Gaussian  beam  with  the  same  waist  width,  the  Af*  factor  can  be  determined.  This 
M1  factor  is  a  fundamental  property  of  the  beam  which  is  not  changed  by  a  lens  if  the  lens 
has  a  sufficient  diameter  to  avoid  truncation  of  the  beam  profile  and  negligible  spherical 
aberration  [8].  In  fact,  the  A/2  determines  the  intensity  at  the  lens  focus  and  the  length 
of  the  focus;  the  product  of  these  quantities  is  a  factor  1/M*  lower  compared  to  a  TEMoo 
beam. 

The  amplification  of  the  Stokes  beam  by  a  nondiffraction-limited  pump  beam  is  solved 
numerically  using  our  numerical  model  described  in  chapter  2.  For  the  nondiffraction- 
limited  pump  beam,  Gaussian- Hermite  (G-H)  beams,  Gaussian-Laguerre  (G-L)  beams,  and 
Gaussian-Scheii-model  (GSM)  beams  are  used.  The  M3  factor  of  these  beams  can  be 
calculated  analytically.  A  random  superposition  of  G-H  modes  is  also  used  as  a  pump 
beam;  for  this  type  of  beam,  the  M*  factor  is  obtained  numerically. 

It  is  shown  that  for  high  amplification  (i.e.  at  threshold  power)  the  gain  of  the  Stokes 
beam  is  approximately  the  same  for  beam  profiles  of  different  type  but  with  the  same  M2. 
Therefore,  the  threshold  power,  i.e.  the  power  at  which  the  Stokes  power  is  1%  of  the  pump 
power,  depends  strongly  on  the  M  2  factor  and  to  a  lesser  extent  on  the  exact  beam  profile  of 
thf  jump.  The  relation  between  M2  and  the  threshold  power  was  confirmed  experimentally 
in  a  Raman  cell  containing  high  pressure  methane. 

Chapter  5  describes  the  optimization  of  the  wavelength  conversion  at  high  pump  energies 
by  the  use  of  astigmatic  lenses  or  a  pair  of  cylindrical  lenses.  Approximately,  the  conversion 
does  not  change  by  varying  the  focal  length  of  an  ordinary  lens  since  the  length  of  the  focus 
times  the  intensity  at  the  focus  is  constant  [9].  Thus,  a  shorter  focal  length  results  in  a 
shorter  Rayleigh  range  (length  of  focus)  which  is  compensated  by  the  increased  intensity 
at  the  focus.  However,  an  astigmatic  focus  can  be  used  to  modify  the  gain  and  therefore 
the  conversion.  An  astigmatic  focus  gives  a  tuning  facility  by  varying  AF,  i.e.  the  distance 
between  the  two  foci. 

In  the  experimental  study,  a  Nd:YAG  laser  was  focused  into  a  Raman  cell  containing  88  bar 
of  methane.  The  conversion  of  the  pump  beam  to  other  wavelengths  was  studied  as  a 
function  of  pump  energy  for  an  ordinary  and  for  an  astigmatic  focus.  We  monitored  the 
energy  at  the  following  generated  wavelengths:  the  first  Stokes  at  1.5  pm,  the  second  Stokes 
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at  2.8  pm,  the  anti-Stokes  at  0.8  /mi,  and  the  Brillouin  scattering  at  1.06  pm.  A  dichroic 
mirror,  transparent  for  the  pump  and  highly  reflective  for  the  first  Stokes,  was  placed  at 
the  entrance  of  the  cell  to  enhance  the  first-Stokes  generation. 

Numerical  simulations  were  performed  and  compared  with  the  experimental  results.  All 
input  parameters  were  obtained  from  the  literature  and  none  of  the  parameters  were  fitted  to 
match  the  simulations  to  the  experiments.  A  reasonable  agreement  between  the  experiment 
and  the  simulation  was  found.  It  is  shown  that  the  increased  conversion  efficiency  for  an 
astigmatic  focus  is  due  to  the  reduction  of  second  Stokes  radiation  at  2.8  pm. 
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2  THEORY  OF  NUMERICAL  SIMULATION 

2.1  Introduction 

A  numerical  model  has  been  developed  at  our  laboratory  to  study  and  to  optimise  wave¬ 
length  conversion  of  laser  radiation  by  stimulated  Raman  scattering  (SRS).  The  extension 
to  an  astigmatic  focus  is  included  in  the  numerical  model  for  the  analysis  of  our  experiments 
with  astigmatic  lenses.  These  experiments  showed  that  the  use  of  an  astigmatic  lens  or  a 
pair  of  cylindrical  lenses  can  increase  the  conversion  efficiency  of  a  Raman  amplifier  [13]. 
In  chapter  5  this  improvement  of  the  conversion  efficiency  for  astigmatic  lenses  is  studied 
further.  In  this  chapter,  we  present  the  numerical  method  in  detail  and  we  show  some 
numerical  results. 

The  Stokes  beam  is  amplified  in  a  spatially  nonuniform  gain  medium  generated  by  the 
focused  pump  beam.  The  problem  of  beam  propagation  in  a  rotationally  symmetric  gain 
medium  has  been  treated  by  expanding  the  pump  and  the  Stokes  beam  in  Gaussian-Laguerre 
modes  and  calculating  the  mode  coupling  [5,6,7].  These  studies  have  shown  that  there  is  a 
narrowing  of  the  Stokes  beam  because  of  the  higher  amplification  at  the  center  of  the  pump 
beam.  This  effect  is  called  gain  focusing  and  the  analytical  model  describing  it  [6]  will  be 
compared  with  our  numerical  results. 

Because  of  the  astigmatic  focus,  the  gain  medium  is  no  longer  rotationally  symmetric. 
In  this  chapter  the  beam  propagation  in  a  non-rotationally  symmetric  gain  medium  is 
studied.  A  method  based  on  finite  differences  is  used  and  this  method  will  be  described  in 
section  2.2.  Since  there  is  no  rotational  symmetry,  both  the  x  and  y  coordinate  have  to 
be  used  instead  of  a  single  radial  coordinate.  Therefore,  a  two-dimensional  grid  is  used  in 
the  finite  difference  (FD)  method.  Because  the  grid  is  propagated  in  the  z-direction,  the 
FD  method  is  in  effect  a  three-dimensional  model. 

The  results  of  the  FD  method  are  cross-checked  using  analytical  models.  It  is  shown  that 
for  low  pump  powers  the  results  of  the  FD  method  agree  with  the  analytical  model  of  Kung 
[14,15].  For  high  powers  the  effect  of  gain  focusing  increases  the  amplification  of  the  Stokes 
beam.  In  the  regime  of  high  pump  powers,  the  results  of  the  FD  method  agree  with  the 
analytical  model  of  Cotter  et  al.  [6]  that  includes  gain  focusing. 

For  the  analysis  of  the  time  dependence  of  the  pump  and  the  Stokes  pulse,  we  have  developed 
a  plane-wave  model  that  uses  one  space  and  one  time  coordinate  without  the  two  transverse 
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space  coordinates.  This  method  uses  the  amplification  obtained  from  the  FD  method,  so  in 
this  way  gain  focusing  and  astigmatism  are  included.  The  details  of  this  method  will  also 
be  presented  in  section  2.2. 

Since  the  experiments  deal  with  the  Stokes  energy  versus  the  pump  energy,  the  experimental 
results  are  compared  with  the  results  of  the  plane  wave  model,  which  include  the  calculated 
energies  of  the  pump  and  the  Stokes  pulse. 


2.2  Beam  propagation  method 

2.2.1  The  paraxial  wave  equation 

The  propagation  of  beams  in  free  space  is  given  by  the  wave  equation 

V*E-  \d*E  =  0,  (2.1) 

c* 

where  E  is  the  scalar  electric  field  and  c  is  the  speed  of  light.  If  we  substitute  in  equation 
(2.1)  the  solution  for  E(x,y,z,t)  of  the  form 

E{x, »,  z,  t)  =  t  ( x ,  v,  z)  expfifcz  -  «'wt]  (2.2) 

and  use  the  slowly  varying  envelope  approximation 

\d\E\  «  \2kdte\,  (2.3) 

then  the  paraxial  wave  equation  is  obtained 

d\C  +  d\e  +  2  ikd.e  =  0 .  (2.4) 

The  relation  u  =  ck  between  the  angular  frequency  w  and  the  wave  number  k  was  used  in 
this  derivation. 

The  paraxial  wave  equation  (2.4)  is  used  for  the  propagation  of  the  Stokes  and  the  pump 
beam.  The  Stokes  field  1,  is  amplified  in  a  gain  medium  that  is  generated  by  the  pump 
field  £p.  The  paraxial  wave  equation  for  the  Stokes  field  is  given  by  [6,16) 

d\e.  +  ale.  +  2  ik,d,e.  =  ik.ge,\eP\* , 


(2.5) 
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where  g  is  the  Raman  gain  coefficient.  The  paraxial  wave  equation  for  the  pump  field, 
including  pump  depletion,  is  given  by 


d\tp  +  d\tv  +  2ik,d.£p  =  -i-*-gtr\lt\' . 


(2.6) 


These  two  partial  differential  equations  can  be  considered  as  an  initial  value  problem  (with 
z  instead  of  t).  The  Stokes  and  the  pump  field  at  the  entrance  window  of  the  Raman 
amplifier  are  known  and  solving  (2.5)  and  (2.6)  gives  us  the  two  fields  at  the  exit  window. 
A  FD  method  is  used  to  solve  (2.5)  and  (2.6)  by  numerically  integrating  these  equations 
along  the  z-axis.  The  numerical  integration  uses  the  values  on  a  two  dimensional  grid  to 
represent  the  cross-section  of  the  Stokes  and  the  pump  field. 

We  shall  introduce  new  transverse  coordinates  because  the  FD  method  needs  a  grid  that 
is  adjusted  to  the  width  of  the  pump  beam.  Otherwise,  only  a  few  grid  points  are  used 
effectively  at  the  focus.  The  new  transverse  coordinates  scale  to  the  width  of  the  pump 
beam  and  are  used  to  define  the  grid.  The  pump  beam  is  focused  into  the  Raman  amplifier 
by  an  astigmatic  lens.  If  there  is  no  pump  depletion  and  the  pump  beam  is  Gaussian,  the 
cross-section  of  the  beam  is  an  ellipse  and  the  widths  of  the  elliptic  Gaussian  pump  beam 
are  given  by  (17} 


w*(z)  =  u£(l  +  (z  -  z/x)1/^.)  "»d 


(2.7) 


W](z)  =  U>*tf(l  +  (Z  -  Zfyf/zly)  , 


(2.8) 


where  in  the  z-z  plane  z/x  is  the  place  of  the  focus,  z0Z  is  the  Rayleigh  range  (2zox  is 
often  called  the  confocal  parameter),  and  woz  is  the  half  width  of  the  beam  at  the  focus. 
The  parameters  for  the  y-z  plane  are  analogous  to  the  parameters  for  the  z-z  plane.  The 
equations  for  the  beam  width  are  used  to  define  two  new  transverse  coordinates  z  and  y 

x  —  x vx(z)x  and  y  =  u/v(z)y.  (2.9) 

The  paraxial  wave  equations  (2.5)  and  (2.6)  can  be  written  in  terms  of  the  new  transverse 
coordinates.  The  equation  for  the  Stokes  field  with  the  new  transverse  coordinates  is 

£.  +  2ik,(dt£,  -  zd.ilnwJdjS. 

*« 

-  ya,(ln  u>y)dgtt)  ~  ik,gt,\£p\2 . 


(2.10) 
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As  a  consequence  of  the  introduction  of  the  new  transverse  coordinates  there  appear  first 
order  transverse  differentials  of  These  differentials  might  be  incorporated  in  a  FD 
method.  However,  we  will  use  a  transformation  that  removes  the  first  order  transverse 
differentials.  The  advantage  of  this  transformation  is  that  the  remaining  term  can  be 
integrated  analytically.  The  first  order  transverse  differentials  of  C,  can  be  removed  by  the 
following  transformation: 


S(x,y,z)  =  £t(x,y,z) 


k,u>l 


exp  [*??^*/*  "  exp  -  *)vJ 

transformed  Stokes  field  S  i 

al  d'±  (  w2  X 2 

~\S  +  —^5  +  ik,  ( da  In  wz  +  d,  In  wv  +  ik, -r  - - — 

*1  \  zl,  +  (*  ~  zlxY 

+  ik»-n  t  — « ) s  +  2ik‘d*s  =  '^yslpl  • 

4, +  (*-*/«)*/ 


(2.11) 


The  equation  for  the  transformed  Stokes  field  S  is 

dl 


(2.12) 


^oyV2 

~oy  1  (Z  —  zIy)i  > 

The  transformation  of  the  pump  field  is  obtained  by  replacing  k,  with  kp  in  the  phase 
factors.  The  equation  for  the  transformed  pump  field  P  is 


wl 


d2  / 

P  +  -~P  +  ikp  ^3,  In  wx  +  ds  In  wy  +  ik 


wLx2 


+  ik 


v*lv  +  (*  -  */v)J 


P*U  +  (z  -  */.)2 

')P  +  2ikpd,P=-i%-gP\S\i. 


(2.13) 


The  equations  (2.12)  and  (2.13)  will  be  solved  numerically  by  the  FD  method  that  is  de¬ 
scribed  below. 


2.2.2  Description  of  the  finite  difference  method 

The  equations  for  the  Stokes  and  the  pump  field  are  solved  by  a  FD  method.  Since  the 
problem  is  an  initial  value  problem,  the  strategy  is  to  integrate  stepwise  along  the  z-axis. 
The  FD  method  will  solve  (2.12)  and  (2.13)  using  the  values  on  a  two  dimensional  grid 
to  represent  the  fields  at  a  cross  section  of  the  beam.  These  grid  values  will  be  used  to 
compute  the  fields  on  these  grid  points  after  some  step  A z  along  the  z-axis.  This  procedure 
starts  at  the  entrance  window  and  ends  at  the  exit  window  of  the  Raman  amplifier. 

The  Stokes  and  the  pump  field  at  the  entrance  window  (z  =  0)  are  represented  by  the  values 
on  a  two  dimensional  grid.  These  grid  values  are  given  by 


«I^°  =  S{x,  y,  z  =  0)  and  p,""0  =  P(z,  y,  z  =  0)  , 


(2.14) 


Figure  2.1:  Rectangular  grid  used  by  the  finite  difference  method.  As  the  grid  advances  along 

the  z-axis  its  width  is  adjusted  to  the  width  of  the  pump  beam. 

with  x  =  l  Ax,  y  =  mAy,  and  z  =  nAz .  (2.15) 


The  subscripts  l  and  m  are  indices  that  indicate  the  position  on  the  grid;  the  superscript 
n  is  an  index  that  gives  the  number  of  integration  steps.  This  grid  is  advanced  along  the 
z-axis  as  is  illustrated  in  figure  2.1.  The  width  of  the  grid  is  adjusted  to  the  width  of  the 
pump  beam  by  using  the  transverse  coordinates  x  and  y  instead  of  x  and  y. 


The  FD  method  for  the  integration  along  the  z-axis  is  composed  of  several  distinct  methods 
obtained  by  splitting  the  differential  equation.  This  procedure  is  based  on  the  FD  operator 
splitting  method  that  is  used  for  initial  value  equations  [18].  The  equation  for  the  Stokes 
field  is  of  the  form  dxS  —  CS,  where  C  is  some  differential  operator.  This  operator  can  be 
split  into  a  linear  sum  of  four  differential  operators 


+ acj** + - 1  (*■*”  ”■ +a-'"' 
+ ,t. ,  ft8*  ..y 

*U  +  (*~  */*)3  *J«  +  (*-*/*)*  / 


—  CiS  +  S  +  CiS  +  C4S . 


(2.16) 


The  FD  method  for  the  equation  dMS  =  CiS  is  called  17,.  With  these  methods  C/,-,  the 
method  for  the  total  operator  is  obtained  by  successively  applying  these  methods  Ui  to  the 
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grid  values.  This  means  that  the  method  for  the  total  operator  is  given  by 

r1  =  I/,(0 
r*  =  U^r1) 
rs  =  Us(ra) 

#"+1  =  I/4  (r8)  (2.17) 

where  the  r*  are  intermediate  grid  values  that  are  discarded  after  completion  of  the  step  in 
z.  The  four  methods  U,  of  (2.17)  will  be  described  below. 


2.2.3  The  method  for  the  diffraction 


The  diffraction  of  the  Stokes  beam  is  given  by  the  second  order  partial  differential  equation 


a,s  =  -^-Ta|5  + 

2k,w\ 


2k.w2 * 


dlS=  thS  +  hS. 


(2.18) 


This  equation  can  be  split  into  two  operators.  The  differentials  in  the  equation  dMS  =  CiS 
are  approximated  by  finite  differences  in  the  following  way 


_n+l  _ 
*i.m 


9l,m 


A  z 


*  y  ik,w2(Ax)2  lktw2(Ax)2  J 


(2.19) 


This  particular  method  of  using  the  present  s'*  and  next  *n+1  grid  values  for  the  second 
order  differential  is  called  the  Crank-Nicholson  scheme  [18,19].  This  scheme  is  stable,  which 
means  that  rapid  fluctuations  of  the  grid  values  that  are  inevitably  introduced  by  round-off 
errors  are  not  amplified  [18].  The  rapid  fluctuations  would  otherwise  suppress  the  important 
slowly  varying  part.  The  following  tridiagonal  matrix  equation  has  to  be  solved  to  obtain 
the  next  grid  values  en+l  from  a" 


~iv9l+l,m  +  (1  +  2»vKm  -  +  (1  "  2 ‘»Km  +  **'*1-1, m .  (2  2°) 

where  v  =  Az/4ktw2{Ax)2.  The  equations  at  the  grid  boundaries  /  =  0  and  l  =  L  are 

-««”?£  +  (1  +  **')*o£  =  *Wi,m  +  (1- and  (2.21) 


ft 

m 


(1  +  iv)«2^m  -  *‘'*2-1.™  =  (1  -  «‘')*2.™  +  *W*2-1,™  • 


(2.22) 
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The  equations  (2.20),  (2.21),  and  (2.22)  can  be  written  in  the  following  matrix  notation 

With  Sn  =  («3,m.«l(m.  • 

(1  4-  »vH)sn+1  =  (1  —  ivH)sn ,  (2.23) 

where  H  is  an  hermitian  matrix  (i.e.  Hy  =  H^).  An  equation  of  the  form  (2.23)  is  unitary 
[20]  meaning  that  |sB+lj  =  |sn|.  This  ensures  conservation  of  power  since  the  power  of  the 
Stokes  beam  is  proportional  to 

The  FD  method  Uj  for  the  solution  of  the  differential  equation  dtS  =  CjS  is  obtained  by 
modifying  method  U\.  A  tridiagonal  matrix  equation  similar  to  (2.23)  can  be  derived  for 
operator  £%,  and  the  method  f/j  consists  of  solving  that  matrix  equation. 

The  pump  field  is  treated  in  a  similar  way  as  the  Stokes  field  except  that  k,  is  replaced  by 

kp. 


2.2.4  The  method  for  the  Raman  amplification 

The  amplification  of  the  Stokes  beam  and  the  depletion  of  the  pump  beam  is  given  by 

dMS  =  |s|P|J  and  (2.24) 

d,P=-  (225) 

These  equations  can  be  solved  in  closed  form  using  the  conservation  of  the  (normalized) 
photon  flux  [21,22]  given  by 

F  =  |S|*  +  «e|P|* ,  (2.26) 


where  k  =  k,/kp.  The  solution  of  (2.24)  and  (2.25)  using  the  constant  photon  flux  (2.26)  is 

I S(z  +  *z)\>  _  \S±z)l  ,  ,  {  , 

\P(t  +  A*)|*  “  |P(*)|*  exPl»FA*/*J  •  l2-27) 


MS 


(2.27) 
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The  FD  method  I/j  for  the  Stokee  and  the  pump  field  is  baaed  on  (2.%)  and  (2.27).  The 
new  values  for  the  Stokes  and  the  pump  field  on  the  grid  are  given  by 


(2.28) 


rt#»+x  _  nn 

Pl,m  ~~  Pi, mi 


X  "ll£Ja  +  KJa«P(*/iVW*]  * 


where 


(2.29) 


(2.80) 


2.2.5 


The  method  for  the  grid  adjustment 


The  introduction  of  the  transverse  coordinates  x  and  y  results  in  the  following  differential 
equation  for  the  Stokes  field: 

dxS  =  —  ^  ^3,  In  wx  +  dg  In  wv 

+  tk,  W;*~ - rr  +  I kt  - rz)  S  .  (2.81) 

This  equation  can  be  solved  in  closed  form.  The  Stokes  field  at  z  +  As  is  given  by 
S(z  +  Az)  =  S(z)i 


.  /  w,(z)wt(z) 
Y  tw,(z  +  Az' 


Az)wv(z  +  Az) 


x  exp 


2  V  zox 


-l  */«  ~  * 
*09 


*oy  *ot  )  u_ 


u=jt+ A# 


(2.82) 


.  he  Stokes  field  is  multiplied  by  an  amplitude  factor  A(z,  Az)  and  a  phase  factor  $(x,  y,  z,  Az). 
The  FD  method  U4  for  the  Stokes  field  is 


*n+l  —  Mn  i»i*n 

vl,m  > 


(2.33) 


where  An  and  are  based  on  (2.32).  The  conservation  of  energy  is  ensured  by  the 
amplitude  factor.  The  power  of  the  Stokes  beam  PJ*  is  proportional  to  £(tm  |z"m|*AxAy. 
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Apart  from  the  Raman  amplification,  the  term  „  is  conserved  for  all  methods  t/». 

The  amplitude  factor  in  the  grid  adjustment  changes  m  |*£m|*  to  compensate  for  the 
change  in  AxAy,  so  that  P?  = 

The  method  for  the  pump  field  is  obtained  by  replacing  k,  with  h,  in  the  phase  factor 

2.2.6  Numerical  stability  and  conservation  of  power 

A  convenient  method  of  testing  the  stability  of  a  FD  method  is  the  von  Neumann  stability 
analysis  [18]  (or  Fourier  analysis  [19]).  If  the  coefficients  of  the  difference  equation  are  con¬ 
stant  (or  very  slowly  varying  and  to  be  considered  as  constants),  the  independent  solutions 
are  of  the  form 

«m  =  £n«tt’mA* ,  (2.34) 

where  Jbx  is  a  wavenumber  of  arbitrary  value  and  £  is  a  complex  number  depending  on  kz. 
This  £  is  called  the  amplification  factor  at  the  wavenumber  kz  since  u^+1  =  £u^. 

Stability  of  the  FD  method  is  determined  from  £  as  a  function  of  kx.  If  there  is  a  particular 
value  of  kx  for  which  |£(fc*)|  >  1,  then  the  FD  method  is  called  unstable.  In  principle  such 
an  unstable  method  might  be  used  for  small  propagation  distances  if  the  initial  function  uj, 
does  not  contain  the  unstable  mode(s).  However,  in  practice  these  unstable  modes  are  high 
frequency  modes  (in  kx)  and  these  modes  are  generated  by  the  algorithm  through  round-off 
errors.  It  then  depends  on  the  accuracy  of  the  number  representation  of  the  computer  how 
many  propagation  steps  can  be  taken  before  the  erroneous  high  frequency  modes,  which 
grow  exponentially  with  |£Jn,  overwhelm  the  stable  low  frequency  inodes.  Thus,  it  is  clear 
that  unstable  methods  should  be  avoided. 

We  will  now  use  the  von  Neumann  analysis  on  the  FD  method  for  the  diffraction,  which  we 
simplify  to 

u,  =  «'u„  (2.85) 

Using  the  Crank- N icholson  scheme,  the  following  FD  method  is  obtained: 

C1  -  <C  _  ■  ( +  «&-!  «mVl  -  1  +  l&t1! 

Ax  \  2(Az)*  2(Ax)* 


(2.86) 
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For  the  amplification  factor  (,  the  following  expression  is  found: 


c  = 


(2.S7) 


It  is  clear  that  j£|  =  1  for  all  values  of  kt  and  for  all  Ax  and  Ax,  so  the  FD  method  (2.36) 
is  unconditionally  stable  meaning  that  the  stability  does  not  depend  on  the  step  size  Ax 
or  the  grid  size  Ax.  Furthermore,  the  method  is  also  unitary  because  |£|  =  1  for  all  kx, 
therefore  power  is  conserved. 

In  the  analysis,  we  did  not  consider  the  influence  of  boundary  conditions  on  stability.  Gen¬ 
erally,  boundary  conditions  do  not  have  a  drastic  influence  on  stability  [16].  In  some  cases, 
boundary  conditions  can  just  tip  the  balance  and  a  matrix  method  should  be  used  [16]. 
We  showed  above  that  the  method  for  the  diffraction  is  unitary  for  the  selected  boundary 
conditions  in  which  case  the  boundaries  operate  as  perfect  mirrors. 

TYansparent  boundaries  are  sometimes  desired,  for  instance  in  our  simulation  of  a  spuri¬ 
ous  window  reflection  in  chapter  4.  In  that  particular  case,  a  divergent  Stokes  beam  was 
amplified  by  a  convergent  pump  beam  and  Stokes  power  had  to  be  removed  at  the  grid 
boundaries.  Intuitively,  it  can  be  understood  that  an  extrapolation  of  the  wavefront  acts 
as  a  transparent  boundary 


um  =  cu"  =>  u"+i  =  cu"  (boundary  condition).  (2.38) 

If  [e|  <  1  and  Imc  >  0,  then  (2.38)  represents  a  transparent  boundary.  This  was  indepen¬ 
dently  discovered  by  Hadley  [20],  who  analyzed  this  method  in  some  detail.  The  algorithm 
is  included  as  an  option  in  the  simulation  program. 

Consistency  is  another  requirement  for  a  FD  method.  A  consistent  (or  compatible)  FD  equa¬ 
tion  converges  to  the  differential  equation  as  Ax,  Az  — ♦  0  [16].  This  is  not  trivial  as  the 
stable  Du  Fort-Frankel  formula  for  the  differential  equation  u*  =  uzc  shows: 

C1  ~  <*m~‘  _  <*m+l  ~  ~  «»’  +  «4-l 

2Az  ~  (Ax)* 

Using  the  Taylor  expansions  around  u£,,  such  as 

«m+1  =  u  +  Az  U,  +  ^y-U„  +  0([Ax]S)  , 


(2.40) 
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and  substituting  the  expansions  in  the  FD  method  (2.39),  we  find 

u,  =  u«,  -  ««  +  0([Ax]*)  +  0([As]*)  •  (2.41) 

If  As  and  Ax  approach  0  with  a  constant  ratio  0  =  As/Ax,  the  Du  Fort-Frankel  formula 
solves  the  differential  equation 

«*  =  u«,  -  0iutM  ,  (2.42) 

which  is  not  equal  to  ua  =  u„!  Thus,  the  FD  method  2.39  solves  a  different  differential 
equation  for  each  ratio  0  =  As/ Ax.  Equation  2.41  also  provides  information  on  the  accu¬ 
racy  of  the  method;  since  Ax  and  As  appear  to  the  second  power  in  the  truncation  terms 
0([Az]2)  and  0([As]2),  the  method  is  called  second  order  accurate  in  x  and  s. 

The  Crank-Nicholson  scheme  for  the  diffraction  equation  is  both  stable  and  consistent.  In 
addition,  it  is  accurate  to  second  order  in  x  and  s.  For  a  linear  initial  value  problem,  the 
following  theorem  holds:  [16] 

Equivalence  theorem  of  Luc:  Given  s  properly  posed  initial  boundary  value  prob¬ 
lem  and  a  finite  difference  approximation  to  it  that  satisfies  the  consistency  con¬ 
dition,  then  stability  is  the  necessary  and  sufficient  condition  for  convergence. 

Thus  the  FD  method  for  the  diffraction  converges  to  the  physical  solution.  However,  there 
is  to  our  knowledge  no  mathematical  theorr  to  that  states  that  our  complete  FD  method 
with  operator  splitting  converges  to  the  physical  solution. 

A  fully-explicit  FD  method  has  been  developed  for  the  diffraction.  The  above  mentioned 
Crank-Nicholson  scheme  is  an  implicit  method,  since  it  uses  the  values  at  the  next  time 
step  u"+l  to  compute  the  differential  utt.  An  explicit  method  uses  un  for  all  differentials 
(except  of  course  u(),  however,  explicit  methods  often  suffer  from  instability.  The  following 
scheme  is  fully  explicit: 

=  j(c. + «;+.) + (i  -  £)< + ^(Ci  -  *c + <-■) 

«  < + + C-Y’> . 


n+J/l 

“w 


•c* 


(*•<») 
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where  c  is  a  small  number  for  stabilisation  of  the  method  and  u"+*^  are  intermediate 
values.  Method  (2.43)  is  stable  for 


e 


1  -  2»*  -  Vl  -  4»* 
2v 


and 


v  = 


2  Az 
(A*)* 


<  0.5. 


(2.44) 


For  small  step  sises,  this  method  is  stable  and  at  least  twice  as  fast  on  the  CONVEX  0230 
as  the  Crank-Nicholson  scheme  (for  the  same  step  sise).  However,  a  serious  limitation  is 
the  small  step  size. 


2.2.7  Accuracy  and  speed 

The  accuracy  of  the  method  was  tested  by  comparing  the  FD-propagation  of  a  Gaussian 
beam  with  the  closed  form  solution  [14].  The  Gaussian  beam  propagated  from  z  =  -L/2 
to  z  =  L/2  with  the  focus  located  at  z  =  0.  We  used  the  symmetry  along  the  x-  and  y-axis 
to  increase  the  speed  by  a  factor  four.  The  grid  consisted  of  40  x  40  points  and  mapped  the 
first  quadrant  (because  of  symmetry). 

The  accuracy  of  the  FD- met  hod  was  estimated  from  the  phase  shift  between  z  =  -L/2  and 
x  =  L/2  which  is  given  in  closed  form  by  [14] 

A^  =  2  tan-1  ~  .  (2-45) 

2zo 

The  FD-method  gave  A4>  =  2.6498  rad,  for  L/2zo  =  4  and  a  step  size  of  Az/2zo  =  0.04; 
this  A^  compared  with  equation  (2.45)  gives  a  relative  error  of  0.07  %.  The  propagating 
beam  suffered  no  power  loss  which  is  an  advantage  over  other  methods  [21]. 

The  speed  of  the  calculation  was  30  msec  per  propagation  step.  The  program  ran  on  a 
CONVEX  C230  and  used  the  vectorization  capacity  of  this  computer.  The  Crank-Nicholson 
method,  which  was  used  for  the  second  order  derivatives,  is  essentially  a  sequential  algo¬ 
rithm.  However,  the  two-dimensional  grid  made  it  possible  to  vectorize  the  d\  differential 
in  the  y-direction  and  vice  versa. 

2.3  Plane-wave  model 

For  the  analysis  of  the  time  dependence  of  the  pump  and  the  Stokes  pulse,  we  have  developed 
a  plane-wave  model  that  uses  one  space  and  one  time  coordinate  without  the  two  transverse 


TNO  rvport 


Page 

22 


•pace  coordinates.  Since  the  experiments  deal  with  the  Stokes  energy  versus  the  pump 
energy,  the  experimental  results  are  compared  with  the  results  of  the  plane  wave  model 
that  include  the  calculated  energies  of  the  pump  and  the  Stokes  pulse. 

The  plane  wave  equations  of  SRS  can  be  found  in  many  textbooks  [1],  so  we  describe  these 
equations  briefly.  In  the  forward  scattering,  the  electric  Stokes  field  £/[*,t)  is  coupled  to 
the  electric  pump  field  £p(*tt)  through  the  Raman  material  vibration  Qf{*,t): 


-dt£f  +  d,£,  =  iK2Q}£, , 

(2.46) 

— dt£P  +  d*£p  —  iKsQf£f  with  Ks  =  -E-jf-Kt , 
vp  v,kt 

(2.47) 

dtQ)  +  TQ)  =  ~iKx£s£; , 

(2.48) 

where  v,  and  vp  are  the  light  speed  of  the  Stokes  and  the  pump  field;  and  T  is  the  Raman 
linewidth  (HWHM).  Dispersion  is  not  taken  into  account  by  the  model,  so  we  will  use  v,  = 
vp  =  v  and  K3  =  kp/k,  Kj.  Backward  scattering  is  described  by  similar  equations,  which 
are  obtained  by  substituting  the  backward  scattered  Stokes  field  £*(z,f)  for  £j,  Q/(z,t)  is 
changed  to  Q»(*,t),  and  the  {1/v  dt+dx)£/  differential  must  be  replaced  by  (l/v3f-3,)£*. 

Steady-state  approximation  of  equations  (2.46)-(2.48)  is  often  possible,  since  SRS  is  a  fast 
process  so  the  material  vibration  Q  is  in  equilibrium  with  the  pump  and  the  Stokes  field. 
Under  steady-state  condition,  we  find  for  the  material  vibration 

<?,  =  =^e,e;  <2 «) 

For  the  pump  and  the  Stokes  field,  the  following  equations  are  obtained  using  (2.49): 

\ditj  +  dx£,  =  ^-£j\£p\2  and  ±dt£p  +  3 x£p  =  -hMl£p\£f\*  .  (2.50) 

These  equations  can  be  simplified  further  by  introducing  the  intensity  of  the  Stokes  (F) 
and  the  pump  (P) 

F  =  |f/|*  and  P  =  \£p\* .  (2.51) 

Note  that  we  use  a  normalized  electric  field  so  that  |f|*  is  equal  to  the  intensity.  The 
equations  (2.50)  now  reduce  to 

^dt F  +  dxF  =  gFP  and  ±dtP+dxP=  -^gFP, 


(2.52) 
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when  $  is  the  Raman 
Stokes  intensity  B,  we 

-dtF  +  d,F  = 

v 

^dtB  +  d.B  = 
-dtp  +  a.p  = 

V 


gain  coefficient  given  by  iKiKj/T.  Adding  the  backward  scattered 


obtain 

gFP, 

(2.53) 

9BP, 

(2.54) 

-!£g{F+B)P. 

(2.55) 

The  method  of  characteristics  is  used  to  solve  the  steady-state  equations  (2.53)-(2.55)  as 
well  as  the  transient  gain  equations  (2.46)-(2.48).  We  will  use  the  following  coordinate 
transformation  (the  characteristics): 


/  =  t>f  +  z  and  b  —  vt  -  2 . 


(2.56) 


It  can  be  shown  that  equations  (2.53)-(2.55)  are  cast  in  the  following  form  using  /  and  6: 


dfF  =  *-FP ,  =  and  dfP  = -^(F  +  B)P .  (2.57) 

2  2  Kg  2 


=  -*£’( 


The  intensity  of  the  forward  Stokes  after  a  propagation  step  A/  is  calculated  by 

nf  +  A/,  b)  =  F(f,  b)eF  with  v  =  c'(F+P)A//1.  (2.58) 

Note  that  this  equation  is  similar  to  the  Raman  amplification  method  in  the  beam  propar 
gation  method.  An  intermediate  pump  intensity,  which  is  given  by 


/*(/,*)  =  />(/,*) 


P(f,b)  +  nf,b) 

P(f,b)+F(f,b)cF’ 


(2.59) 


is  used  for  the  calculation  of  the  backward  scattered  Stokes: 

B{f,b+  A6)  =  B(f,  6)cb  with  .  (2.60) 

The  pump  at  (/  +  A/,  6)  is  given  by 


P(f  +  A/,6)  _  — 


(2.61) 
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Equations  (2.58)-(2.61)  can  be  described  by  two  methods,  one  for  the  pump  ( Mf )  and  the 
other  for  the  scattered  radiation  (M*),  that  are  applied  sequentially: 

F(/ +  A/,6)  =  M,(F(/, 6),  P(/, 6))  and  P'{f,b)  =  Mp(F{f,b),P{f,b))  (2.62) 

and  subsequently 

B(/,  6 +  A6)  =  M,(B(/,  6), /*(/,*))  and  P(/  +  A/,6)  =  M„(B(/,  6), /*(/,*)) 

(2.63) 

By  adding  similar  equations  to  this  chain,  other  processes  such  as  second  Stokes  generation 
can  be  added  to  the  simulation. 

Since  the  Raman  medium  is  filled  with  the  pump  radiation  and  the  scattered  radiation  over 
the  total  length  of  the  medium,  the  radiation  is  represented  by  arrays  that  give  the  intensity 
at  a  position  z  —  mAz  at  a  time  t.  The  pump  array  Pm(t)  is  related  to  P[f,b)  by 

Pm(t)  =  P(/,  b)  =  P(vt  +  mAz,  vt  -  mAz) .  (2.64) 

If  z  is  advanced  by  A z,  then  t  should  be  advanced  by  At  =  Az/v  in  order  to  keep  the 
backward  coordinate  b  constant.  Under  this  condition,  P(f  +  A/,6)  is  given  by 

P(f  +  A/,6)  =  Pm+i(t  +  At)  with  A  /  =  2wAt  =  2Az .  (2.65) 

The  amplification  chain  (2.62)-(2.63)  can  now  be  written  in  terms  of  the  intensity  arrays: 

Fm+1(t  +  At)  =  Af,(Fm(t),Pro(t))  and  /*(t)  =  Afp(Fm(t),Pm(t)) ,  (2.66) 

and  subsequently  using  At  =  A/ 

Bm-i(t  +  At)  =  M,(Bm(t),  O))  P>m(t)) .  (2.67) 

In  the  transient-gain  situation,  the  notation  gets  more  complicated  and  the  material  vibra¬ 
tion  has  to  be  included  in  the  simulation.  However,  the  principles  for  the  transient-gain 
simulation  are  identical  to  the  steady-state  situation  described  above. 

We  have  to  modify  the  plane-wave  approximation  to  the  realistic  situation  of  focused  beams. 
In  the  equations  above,  the  intensity  was  used,  which  can  be  transformed  to  the  power 
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divided  by  the  beam  area  A.  The  intensity  equations  (2.53)-(2.55)  can  now  be  written  in 
terms  of  powers  and  beam  area: 

-dtF  +  d,F  =  =  AfP,  (2.68) 


-dtF  +  d,F  =  =  AFP, 

v  A 

-dtB  +  d,B  =  *rBP  =  ABP, 
v  A 

-dtP  +  dtP  =  -%-^[F  +  B)P  =  -!£a(F+B)P. 

V  Kg  Kg 


(2.69) 

(2.70) 


using  the  symbols  F,  B,  and  P  for  the  power  instead  of  the  intensity.  Focused  beams  can 
now  be  represented  by  a  beam  area  A(z)  as  a  function  of  s,  where  A(z)  is  given  by  the  area 
of  the  undepleted  pump.  The  cross-sectional  area  of  the  pump  beam  is 


A(z)  =  xws{z)wv(z) , 


(2.71) 


where  2wz  and  2u>„  are  the  beam  diameters  in  the  x-  and  y-direction,  referring  to  the  1/e2 
intensity  points. 

A  higher  accuracy  is  obtained  by  using  the  amplification  obtained  from  the  beam  propaga¬ 
tion  method  described  in  section  2.2.  The  amplification  of  the  plane  wave  model  is  (without 
pump  depletion) 


expf  f  A(z)Pdz\. 

'J  entrance  ' 


(2.72) 


By  comparing  this  amplification  with  the  amplification  of  the  beam  propagation,  which  is 
given  by  (Fexit/ P«ntnmc«)>  we  can  define  a  gain  focusing  factor  q  by 


_  ln(P«xit/ intranet) 

SSL*  *A{z)Pdz 


(2.73) 


Gain  focusing  is  the  narrowing  of  the  Stokes  beam  due  to  the  higher  amplification  at  the 
center  of  the  pump  beam.  Typical  values  for  rj  are  between  1  and  2.  The  modified  gain 
coefficient  A(z)  is 


A(z)  = 


A{z)- 


(2.74) 


Spontaneous  noise  has  to  be  incorporated  in  the  model  for  the  onset  of  the  scattered  ra¬ 
diation.  In  the  equation  of  the  forward  Stokes,  the  spontaneous  power  6  is  included  as 
follows 


ia,F’  +  3,F  =  A(/,  +  fi)P. 


(2.75) 


TNO  report 


Page 

26 


In  the  simulation,  the  spontaneous  power  is  simply  added  to  equations  (2.66)  and  (2.67): 


Fm+l(t  +  At)  =  M.{Fm{t),  />m(t))  +  A(z)8Pm(t)&z 


(2.76) 


and 


Bm_i(t  +  At)  =  /*,(*))  +  A(*)*Pm(t)Az .  (2.77) 

For  the  transient  gain  situation,  spontaneous  noise  is  introduced  in  the  material  vibration, 
thus  spontaneous  scattering  is  generated  indirectly  through  the  material  vibration. 

Reflections  at  the  windows  of  the  Raman  cell  have  a  significant  influence  on  the  SRS  process. 
Because  of  the  very  high  gain  (e3S  at  threshold),  a  small  reflection  can  still  decrease  the 
threshold  energy  [25].  Furthermore,  a  dichroic  mirror  is  used  to  reflect  the  backward  Stokes 
into  the  Raman  cell  so  that  all  Stokes  radiation  leaves  the  cell  in  the  forward  direction.  In 
the  simulation,  the  reflections  can  be  easily  incorporated  using 

F0(t  +  At)  =  (2.78) 

at  the  entrance  window  and 

BM(t  +  At)  =  R^tFM(t)  (2.79) 

at  the  exit  window.  There  is  no  provision  for  pump  reflections  in  the  model. 


2.4  Results  and  Discussion 

2.4.1  Amplification  in  a  spherical  lens  focus 

The  FD  method  that  was  described  above  is  used  to  obtain  the  amplification  of  a  Gaussian 
Stokes  beam  by  a  Gaussian  pump  beam.  The  pump  beam  is  focused  by  an  ordinary  spherical 
lens  at  the  middle  of  the  Raman  cell  of  length  L.  For  the  Stokes  beam  at  the  entrance 
window,  a  width  and  a  divergence  that  are  identical  to  those  of  the  pump  beam  are  used. 
The  Stokes  wave  number  is  identical  to  the  pump  wave  number  in  the  following  calculations. 
This  particular  problem  has  been  analysed  by  analytical  methods  and  the  results  of  these 
methods  are  compared  with  the  results  of  the  FD  method. 


Table  2.1:  Amplification  of  a  Gaussian  Stokes  beam  by  a  Gaussian  pump  beam  with  L/2z0  = 

4. 


pump  power 

FD  method 

ln[A(£)/P,(0)| 
Hung’s  method  parabolic 
\P p  tan-1  pump 

matched  mode 

{Pp-2y/pp)  tan-1^ 

10 

7.2 

6.6 

6.6 

4.9 

20 

16.3 

13.3 

15.8 

14.7 

40 

37.3 

26.5 

36.9 

36.3 

60 

59.9 

39.8 

59.4 

59.0 

The  following  expression  for  the  amplification  of  the  Stokes  beam  has  been  given  by  Kung 
[14,15] 


In 


P.{L) 
P.(  o) 


^  tan-1 
2 


(2.80) 


where  Pt  and  Pp  are  the  normalized  Stokes  and  the  normalized  pump  power,  and  za  is  the 
Rayleigh  range  of  the  pump.  Depletion  of  the  pump  beam  is  not  considered  in  (2.80).  The 
following  normalization  for  the  Stokes  and  the  pump  power  Is  used  [6] 


P  =  9\£o\2kwl, 


(2.81) 


so  that  P  is  a  dimensionless  quantity;  i0  is  the  electric  field  at  the  focus  in  the  center  of 
the  beam,  and  w0  is  the  half  width  of  the  beam  at  the  focus.  The  amplification  given  by 
(2.80)  for  several  pump  powers  is  listed  in  table  2.1  together  with  the  results  of  the  FD 
method.  Comparing  the  amplifications,  we  find  a  higher  amplification  for  the  FD  method 
and  this  difference  increases  with  higher  pump  power.  This  higher  amplification  for  the 
FD  method  can  be  explained  by  considering  the  narrowing  of  the  Stokes  beam  due  to  the 
higher  intensity  at  the  center  of  the  pump  beam.  The  narrowing  at  high  pump  powers  is 
called  gain  focusing  [5,6,7]. 

We  shall  analyze  the  gain  focusing  effect  using  the  approximation  that  the  transverse  profile 
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of  the  pump  beam  is  represented  by  a  parabola.  The  Gaussian  pump  beam  is  given  by 


wa 


I*’1  = 


**  +  »> 

'  «4(«)  I  • 


(2.82) 


where  wp(z)  is  the  half  width  of  the  beam.  Using  the  parabolic  approximation,  the  gain 
profile  is  given  by 


The  Stokes  beam  is  given  by  [6,17] 

£t  =  exp  ^(*)  +  —  g  y  <?(z))]  ’ 


(2.83) 


(2.84) 


where  <f>(z)  and  Q(z)  are  complex  functions  representing  the  amplitude  and  the  width  of 
the  Stokes  beam.  The  expression  for  the  Stokes  beam  and  the  gain  profile  are  substituted 
in  (2.5),  the  paraxial  wave  equation,  and  the  terms  with  the  factor  x1  +  y*  are  separated 
from  the  terms  without  this  factor,  which  leads  to  the  following  set  of  ordinary  differential 
equations: 


dA 


and 


(2.85) 


(2.86) 


These  two  equations  are  solved  by  the  Runge-Kutta  algorithm. 

The  amplification  given  by  the  parabolic  pump  model  is  listed  in  table  2.1.  It  is  clear 
from  table  2.1  that  the  values  for  high  pump  power  correspond  very  well  with  the  FD 
method.  This  can  be  understood  by  realising  that  the  width  of  the  Stokes  beam  is  narrow 
for  high  pump  power  due  to  the  gain  narrowing.  For  a  narrow  Stokes  beam  the  parabolic 
approximation  is  a  better  approximation  than  for  a  wide  Stokes  beam,  and  this  leads  to 
a  better  correspondence  with  the  FD  method  for  high  pump  powers.  The  effect  of  gain 
narrowing  on  the  width  of  the  Stokes  beam  is  shown  in  figure  2.2.  The  amplification  for  low 
power  agrees  with  the  amplification  given  by  (2.80),  since  for  low  power  the  gain  narrowing 
effect  is  negligible. 
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z-oxis  (z/2z0) 


Figure  2.2:  Width  of  the  Stokes  beam  along  the  *- axis  for  several  pump  powers. 

An  analytical  model  for  the  amplification  including  gain  focusing  has  been  derived  by  Cotter 
et  al.  [6]  using  the  gain  focusing  formulation  of  Kogelnik  [5]  and  calculating  the  amplifi¬ 
cation  of  the  “matched”  Gaussian  mode.  The  matched  mode  is  the  mode  for  which  the 
beam  parameters  remain  constant  while  propagating  through  the  gain  medium  [5,6,7].  The 
amplification  of  the  matched  mode  is  given  by  [6] 

=  <2'87» 

In  table  2.1  the  results  of  (2.87)  for  several  pump  powers  are  shown.  The  results  of  (2.87) 
correspond  very  well  with  the  FD  method  for  high  pump  powers. 


2.4.2  Amplification  in  an  astigmatic  focus 

We  shall  use  the  FD  method  to  obtain  the  amplification  of  an  elliptic  Gaussian  Stokes  beam 
by  an  elliptic  Gaussian  pump  beam.  The  pump  beam  is  focused  by  an  astigmatic  lens  or 
by  a  pair  of  cylindrical  lenses  at  the  middle  of  the  Raman  cell.  The  astigmatic  focus  of 
a  Gaussian  beam  is  characterized  by  (2.7)  and  (2.8),  where  z/x  and  z/v  are  the  place  of 
the  foci  in  the  x-z  and  the  y-z  plane,  respectively.  An  important  parameter  describing  the 
astigmatic  focus  is  the  distance  A F  between  the  two  foci,  which  is  | z/,  -  z/J.  Figure  2.3 


t 


m 


■vm* 
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z-oxis  (2/2 z0) 


Figure  2.3:  Intensity  on  the  optical  axis  for  several  values  of  A F. 

shows  for  several  values  of  A F  the  intensity  of  the  pump  beam  on  the  optical  axis  given  by 

m 


IPix  —  °» !/  =  0,  z)  =  /0 


Vlx(z)Wy(z)  ’ 


(2.88) 


where  I0  is  a  constant;  wz  and  tvv  are  given  by  (2.7)  and  (2.8).  It  is  clear  that  the  intensity 
profile  is  strongly  influenced  by  A  F. 

The  amplification  of  the  Stokes  beam  versus  A F  is  shown  in  figure  2.4.  Obviously,  the 
amplification  decreases  for  increasing  A F.  This  decrease  of  the  amplification  is  caused  by 
a  decrease  of  the  intensity  in  the  focal  region.  In  figure  2.5  the  intensity  at  the  focus  in  the 
center  of  the  beam  versus  A F  is  shown,  and  it  is  clear  that  the  intensity  drops  rapidly  as 
A F  is  increased. 

At  high  light  intensity,  Raman  media  are  transformed  into  a  high  absorbing  plasma.  This 
transition  into  the  opaque  state  is  called  optical  breakdown  [26,27].  This  breakdown  effect 
has  to  be  suppressed  in  an  efficient  Raman  cell.  Therefore,  the  power  of  the  pump  beam 
is  limited  so  that  the  intensity  at  the  focus  does  not  exceed  the  breakdown  threshold. 
However,  decreasing  the  power  of  the  pump  beam  will  also  reduce  the  amplification  of  the 
Stokes  beam. 
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Figure  2.4:  Amplification  of  the  Stoke*  beam  ver*u.  A F  tor  several  pump  power*.  The  fol 

lowing  parameter*  are  used:  *o  ~  *o%  *  *o»>  *  *p»  4n“  L/2x0  —  10 ■ 


>■  0.4 
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*F/2*. 


Figure  2.5:  Intensity  at  the  focu*  in  the  center  of  the  pump  beam  versus  AF. 
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Astigmatic  lonaaa  can  be  used  to  reduce  the  intensity  in  the  focus  and  thus  eliminating 
optical  breakdown  without  increasing  the  length  of  the  Raman  cell.  As  figure  2.5  shows, 
the  intensity  in  the  focus  drops  rapidly  with  increasing  AF;  a  AF  equal  to  the  confocal 
parameter  2m,  reduces  the  intensity  with  a  factor  2.  Figure  2.4  shows  that  the  amplification 
is  reduced  only  slightly  with  a  factor  1.2  for  AF  —  2m,.  Thus,  an  astigmatic  lens  gives 
approximately  the  same  amplification  as  a  spherical  lens,  but  with  a  lower  intensity  at  the 
focus. 

For  a  high  energy  pump  beam  in  combination  with  a  tight  focus,  the  conversion  to  a  high 
energy  Stokes  beam  is  limited  by  the  number  of  atoms  or  molecules  in  the  ground  state 
[1,22].  This  saturation  effect  is  called  ground-state  depletion  and  occurs  when  the  relaxation 
time  of  the  final  Raman  level  is  long  compared  to  the  pulse  width  of  the  laser.  Under  these 
conditions  the  number  of  generated  Stokes  photons  is  limited  by  the  number  of  atoms  (or 
molecules)  in  the  volume  where  the  conversion  takes  place,  i.e.  the  focal  region.  We  calculate 
the  volume  of  the  focal  region  using  the  area  of  the  cross-section  of  the  elliptic  Gaussian 
beam  and  the  length  of  the  focal  region: 

V  =  f  xwx{z)wy(z)  dz ,  (2.89) 

where  ws  and  wt  are  the  widths  of  the  beam  given  by  (2.7)  and  (2.8);  Zi  and  zt  are  the 
positions  on  the  optical  axis  where  the  intensity  is  half  the  maximum  intensity.  Figure  2.6 
shows  the  volume  given  by  (2.89)  as  a  function  of  AF.  It  is  clear  that  the  volume  increases 
rapidly  with  A F;  for  AF  =  2 z0  the  volume  is  about  a  factor  3  larger  compared  to  AF  =  0. 
This  means  that  the  energy  of  a  Stokes  beam  limited  by  ground-state  depletion  can  increase 
substantially  by  using  an  astigmatic  focus.  However,  chapter  5  shows  that  ground-state 
depletion  is  negligible  in  our  experimental  setup.  Furthermore,  simulations  show  that  the 
effective  volume  is  significantly  larger  than  the  volume  given  by  (2.89). 

In  a  Raman  cell  with  spherical  lenses,  the  intensity  at  the  focus  can  be  reduced  and  the 
volume  of  the  focal  region  can  be  increased  by  the  use  of  a  longer  focal  length.  However,  this 
approach  has  the  disadvantage  that  it  results  in  an  unwieldy  long  Raman  cell.  Furthermore, 
the  length  of  a  Raman  cell  is  not  as  easily  adjusted  as  the  distance  between  two  cylindrical 
lenses  (changing  AF).  Another  advantage  of  a  short  Raman  cell  with  an  astigmatic  focus  is 
related  to  the  gain  enhancement  by  a  dichroic  mirror  at  the  entrance  window.  The  function 
of  such  a  dichroic  mirror  is  to  reflect  the  backward  scattered  Stokes  beam  back  into  the  cell 
while  not  affecting  the  pump  beam  [14,28].  A  short  Raman  cell  is  essential  for  this  method 
of  gain  enhancement  since  the  delay  of  the  reflected  Stokes  pulse  must  be  short  in  order  to 
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Figure  2.6:  Volume  of  the  focal  region  versus  AF. 


overlap  with  the  pump  beam. 


2.5  Summary  and  Conclusions 

A  three-dimensional  finite  difference  method  has  been  presented  for  the  propagation  of 
beams  focused  by  an  astigmatic  lens.  This  FD  method  is  stable  and  unitary,  which  means 
that  the  power  of  the  beams  is  conserved. 

The  FD  method  is  used  for  the  modeling  of  a  Raman  amplifier.  The  results  of  the  FD 
method  have  been  cross-checked  with  analytical  models  and  with  a  numerical  model  that 
describes  gain  focusing.  It  has  been  shown  that  the  FD  method  gives  reliable  results  both 
for  low  and  for  high  pump  power. 

At  high  light  intensity,  optical  breakdown  occurs  and  this  effect  can  limit  the  efficiency 
of  a  Raman  cell.  It  has  been  shown  that  astigmatic  lenses  lessen  the  problem  of  optical 
breakdown  by  reducing  the  intensity  at  the  focus  while  maintaining  sufficient  gain. 
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3  EXPERIMENTS  USING  METHANE 

3.1  Experimental  setup 

A  diagram  of  the  experimental  setup  is  shown  in  fig.  3.1.  A  Q- twitched  Nd.YAG  laser  was 
used  to  generate  the  pump  beam.  The  BM1-501-D.NS  laser  provided  1.06  ftm  pulses  up  to 
350  mJ  in  energy  with  a  F  WHM  pulse  width  of  7  ns.  Using  a  fast  photo  detector  with  a  rise 
and  fall  time  of  1  ns,  we  measured  the  pump  power.  The  pump  pulses  showed  some  beating 
as  can  be  seen  from  the  1  GHz  oscilloscope  trace  shown  in  fig.  3.2.  Beating  at  300  MH*  can 
be  attributed  to  the  longitudinal  mode  spacing  of  the  50  cm  long  laser  resonator.  A  circular 
diaphragm  could  be  placed  inside  the  laser  cavity  to  control  the  quality  of  the  beam.  We 
used  two  diaphragm  sixes  in  our  experiments,  namely  a  diameter  of  1.5  mm  and  a  diameter 
of  6.5  mm.  The  beam  quality,  i.e.  the  M*  factor,  was  measured  using  a  standard  knife-edge 
technique;  details  of  the  M1  measurement  are  described  in  section  3.3.  The  energy  of  the 
laser  beam  was  adjusted  by  an  attenuator  before  the  beam  entered  the  Raman  cell. 

The  pump  beam  was  focused  into  the  Raman  cell  by  an  ordinary  spherical  lens  of  10  cm 
focal  length  or  by  a  pair  of  cylindrical  lenses  of  10  cm  focal  length  placed  1.3  cm  apart. 
At  the  entrance  of  the  Raman  cell  a  dichroic  mirror  was  placed  to  reflect  the  backward 
scattered  Stokes  into  the  Raman  cell  in  order  to  obtain  all  Stokes  radiation  in  the  forward 
direction.  For  the  threshold  measurements  of  chapter  4,  we  did  not  use  the  dichroic  mirror 
and  we  only  used  the  spherical  lens. 

The  20  cm  long  Raman  cell  contained  methane  with  a  pressure  of  88  bar.  A  window  was 
mounted  on  top  of  the  cell  so  that  visible  sparking  due  to,  unwanted,  optical  breakdown 
could  be  observed.  To  avoid  optical  breakdown,  we  measured  upto  200  mJ  for  the  ordinary 
focus;  for  the  astigmatic  focus  we  measured  upto  250  mJ.  No  sparking  was  observed  below 
these  energies. 

Calibrated  pyroelectric  detectors  were  used  to  measure  the  pulse  energies  of  the ,  amp  pulse 
and  the  generated  pulses.  Various  filters  at  the  exit  window  (selective  mirrors,  color  glass 
filters)  were  used  to  select  and  measure  the  first  Stokes,  the  second  Stokes,  or  the  anti- 
Stokes.  The  energy  of  the  initial  pump  beam  was  measured  with  a  second  detector  using 
the  reflection  of  the  beam  sampler.  Brillouin  scattering  was  measured  by  monitoring  the 
backward  scattered  1.06  (im  beam. 
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8.2  Methane  parameters 


For  our  simulations,  the  Raman  gain  coefficient  has  to  be  obtained.  The  expression  for  the 
gain  of  SRS  is  given  by  (26] 


2X*AN  (do\ 

*  hv,*A v  \5n/  ’ 


(8.1) 


where  A,  is  the  Stokes  wavelength  in  the  medium;  hv,  is  the  energy  of  the  Stokes  photon; 
AN  is  the  population  difference  between  the  initial  and  final  states;  Av  is  the  linewidth 
(FWHM);  and  (do/dZl)  is  the  differential  scattering  cross  section  in  intensity  ratio.  Since 
these  quantities  depend  on  the  pump  wavelength  and  on  the  gas  pressure,  the  Raman  gain 
coefficient  is  calculated  using  the  following  relations: 


•  {do/d fl)  as  a  function  of  pump  wavelength. 

•  dependence  of  Av  on  gas  pressure. 

•  refractive  index  n,  versus  gas  pressure. 

•  molecular  density  Ntot  as  a  function  of  gas  pressure. 


The  differential  scattering  cross  section  is  described  by  [27] 


(3.2) 


where  (do/dfl)^ji9  is  the  differential  scattering  cross  section  of  nitrogen  at  v£;  vp  and  v, 
stand  for  the  pump  and  the  Stokes  frequency;  C(vp,i/p)  is  a  correction  factor  to  account 
for  the  difference  in  the  tabulated  and  the  required  pump  frequency;  and  S(i^)  is  the 
tabulated  relative  cross  section  given  by  Schrotter  and  Klockner  [28].  Since  there  is  no  clear 
wavelength  dependence  of  E(j/')  for  methane,  we  will  use  the  average  value  8.56  for  £  (and 
no  correction  factor).  At  the  pump  wavelength  of  1.06  fan,  a  (do/d fl)  of  7.6x  10~sl  em’/sr 
is  obtained  using  [  (da/dflJ^W)-4]  =  5.05  x  10-48  cm ®/sr,  which  is  constant  for  all  i/t. 
The  accuracy  of  (do/dtl)  is  around  10%,  which  makes  it  the  quantity  with  the  lowest 
accuracy  [26]. 

According  to  Ottusch  and  Rockwell  [26],  the  linewidth  as  a  function  of  the  gas  pressure  is 


Av  =  8.4  x  10®  +  3.9  x  108  p  [Hi] 


with  p  in  [atm]. 


(8.3) 
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So  for  88  bar,  we  obtain  a  linewidth  of  42.7  GHs. 

For  the  refractive  index  of  methane  venue  the  gas  pressure,  we  find  [29} 

n,  =  1  +  4.44  x  10 ~*p  with  p  in  [atm].  (3.4) 

In  the  calculation  of  the  molecular  density  of  methane,  nonideality  of  the  gas  is  accounted 
for  by  using  the  following  expression: 

( P+~){V-nb)  =  nBT ,  (3.5) 

where  n  is  the  number  of  moles,  R  is  the  gas  constant,  and  T  is  the  absolute  temperature. 
For  methane,  the  constants  a  and  6  are  [29] 

a  =  2.253  x  106  cmflatm/mol*  and  b  =  42.78  cms/mol .  (3.6) 

At  a  pressure  of  88  bar,  a  density  Ntot  of  X  53  x  10J1  cm-*  is  found.  For  methane,  we  have 
AN  =  Ntot  (cf.  AN  =  0.66tftot  for  hydrogen). 

Using  these  values,  we  obtain  a  gain  coefficient  of  0.49  cm/GW  for  a  pump  wavelength  of 
1.06  pm  and  a  gain  coefficient  of  0.27  cm/GW  for  the  second  Stokes  generation.  Equar 
tion  (3.1)  shows  that  the  gain  coefficient  is  proportional  to  l/i/„  since  (da/dCl)  is  propor¬ 
tional  to  i/*.  The  relation  between  the  gain  and  the  gas  pressure  is  depicted  in  figure  3.3 
showing  a  saturation  at  high  pressures  due  to  the  increase  of  the  linewidth. 

An  estimation  of  the  spontaneous  Raman  noise  is  made  using  the  spontaneous  intensity 
rhu),/4A  as  derived  in  [30].  A  noise  power  6  of  10-9  W  is  obtained  in  this  way.  Correcting 
for  the  higher  divergence  of  the  pump  beam  with  M1  =  8.7,  we  obtain  a  noise  power  6  of 
M*  x  10-9  W,  i.e.  10"7  W. 


3.3  Beam  quality  measurements 

The  beam  quality  (i.e.  the  M7  factor)  was  measured  using  a  standard  knife-edge  technique. 
A  lens  with  a  75  cm  focal  length  was  placed  in  the  pump  beam,  and  the  width  of  the  beam 
at  various  distances  from  the  lens  was  measured  with  the  knife-edge  technique.  We  used 
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Figure  3.3:  Raman  gain  coefficient  of  methane  versos  gas  pressure  for  a  pump  wavelength  of 

1.06  fim. 

the  clip  levels  of  10%  and  90%,  as  discussed  by  Siegman  in  ref.  [31],  to  obtain  the  beam 
radius  w(z)  using 


w  =  0.78 Dc , 


(3.7) 


where  De  is  the  clip  width,  i.e.  the  separation  of  the  10%  and  the  90%  cut-off  levels.  We 
also  used  the  clip  levels  of  16%  and  84%  using  the  following  expression: 


w  =  1.0 D[ ,  (3.8) 

with  as  the  new  clip  width.  The  difference  for  these  two  methods  was  typically  5%.  Fig¬ 
ure  3.4  shows  the  measured  beam  radius  versus  the  knife  position  as  well  as  the  theoretical 
relation 


»’(«)  =  <4(1  +  ^r^)  with  =  (»•») 

This  expression  for  the  Rayleigh  range  zjt  was  used  to  obtain  the  M*  factor.  For  the  1.5  mm 
diaphragm,  an  M*  factor  of  1.12  ±  0.05  was  found.  Since  this  Ms  value  is  close  to  one,  the 
beam  for  the  1.5  mm  diaphragm  is  almost  diffraction-limited.  A  much  higher  Af2  value  was 
found  for  the  6.5  mm  diaphragm,  namely  8.7  ±  0.3. 


M 
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Figure  3.4:  Beam  radius  versus  distance  from  focus.  The  lines  indicate  the  theoretical  relation 

for  M7  =  1.12  and  M*  =  8.7. 

A  pyroelectric  matrix  array  of  32  x  32  elements  was  used  to  measure  the  pump  beam  profile. 
Figure  3.5  shows  the  near-field  profile  of  the  pump  beam  for  the  6.5  mm  diaphragm.  It  is 
clear  that  the  profile  is  approximately  Gaussian;  i.e.  there  are  no  sharp  peaks  due  to  higher 
order  modes. 


U 
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4  THRESHOLD  OF  STIMULATED  RAMAN  SCATTERING 

4.1  Introduction 

In  this  chapter  we  study  the  effect  of  pump  beam  quality  on  the  threshold  power.  The  SRS 
process  has  no  real  threshold  like  a  laser,  but  an  experimental  threshold  can  be  defined.  Here 
the  threshold  stands  for  the  pump  power  corresponding  to  1  percent  conversion  efficiency. 
By  changing  the  diaphragm  inside  the  laser,  we  could  change  the  quality  of  the  output  beam. 
Because  the  original  laser  beam  (called  the  pump  beam)  is  focused  into  the  Raman  medium, 
the  SRS  process  will  depend  on  the  beam  quality  since  the  process  is  driven  by  the  intensity 
at  the  focus.  We  present  theoretical  and  experimental  results  for  a  nondiffraction-limited 
pump  beam. 

The  beam  quality  is  related  to  the  M 2  factor  which  is  a  dimensionless  factor  that  compares 
the  divergence  of  a  given  beam  with  the  divergence  of  a  (diffraction-limited)  Gaussian  beam. 
The  divergence  of  an  arbitrary  beam  can  be  defined  unambiguously  using  the  variances  a2 
and  o*  of  the  beam  profile  in  the  two  transverse  di.  ections.  It  has  been  shown  that  these 
variances  a\  and  o*  increase  proportionally  to  the  square  of  the  distance  from  the  waist, 
i.e.  where  the  width  is  smallest,  which  is  similar  to  the  width  change  of  Gaussian  beams 
[32,31].  By  comparing  the  rate  of  increase  of  the  variance  of  an  arbitrary  beam  (i.e.  the 
divergence)  with  the  corresponding  rate  of  a  Gaussian  beam  with  the  same  waist  width, 
the  Ml  factor  can  be  determined.  This  M2  factor  is  a  fundamental  property  of  the  beam 
which  is  not  changed  by  a  lens  if  the  lens  has  a  sufficient  diameter  to  avoid  truncation  of 
the  beam  profile  and  negligible  spherical  aberration  [8].  In  fact,  the  M 2  determines  the 
intensity  at  the  lens  focus  and  the  length  of  the  focus;  the  product  of  these  quantities  is  a 
factor  1/M2  lower  compared  to  a  TEMqo  beam. 

The  amplification  of  a  Raman  shifted  beam  (i.e.  the  Stokes  beam)  by  a  nondiffraction- 
limited  pump  beam  is  solved  numerically.  A  finite  difference  (FD)  method  is  used  for  this 
problem  of  beam  propagation  in  a  gain  medium  that  is  generated  by  the  pump  beam. 
The  details  of  the  FD  method  have  been  reported  in  ref.  [33].  For  the  nondiffraction- 
limited  pump  beam,  Gaussian-Hermite  (G-H)  beams,  Gaussian-Laguerre  (G-L)  beams,  and 
Gaussian-Schell-model  (GSM)  beams  are  used.  The  M2  factor  of  these  beams  can  be 
calculated  analytically.  A  random  superposition  of  G-H  modes  is  also  used  as  a  pump 
beam;  for  this  type  of  beam,  the  M2  factor  is  obtained  numerically. 


TNO  report 


Experiments  were  performed  to  confirm  the  relation  between  threshold  and  beam  quality. 
A  NdrYAG  laser  beam  was  focused  into  a  Raman  cell  containing  high  pressure  methane  as 
the  Raman  medium.  The  pump  beam  quality  was  varied  by  changing  the  diaphragm  inside 
the  Nd:YAG  laser  cavity.  Subsequently,  the  Stokes  energy  was  measured  as  a  function  of 
pump  energy.  We  compared  the  experimental  Stokes  energy  with  the  simulated  Stokes 
energy  for  low  pump  energies.  A  good  agreement  was  found  both  for  high  (M1  =  1.12)  and 
low  (Af*  =  8.7)  beam  quality. 


4.2  Theory 


4.2.1  The  beam  propagation  method 


In  the  paraxial  approximation  the  propagation  of  the  Stokes  field  £,  and  the  pump  field  £p 
is  given  by  [6,13] 

d\£,  +  d\£,  +  2  ik,d,£,  =  iktg£, \£,\*  and  (4.1) 


d\£p  +  d%  +  2 ikpdx£p  =  -i%£p\£t\*  , 


(4.2) 


where  k,  and  kp  are  the  wave  numbers  of  the  Stokes  and  the  pump  field  and  g  is  the  Raman 
gain  coefficient.  These  two  partial  differential  equations  can  be  considered  as  an  initial 
value  problem  with  the  space  coordinate  z  instead  of  the  time  coordinate.  The  spontaneous 
Stokes  is  incorporated  by  using  an  initial  seed  beam  which  is  amplified  by  the  gain  medium 
generated  by  the  pump  beam.  A  finite  difference  (FD)  method  is  used  to  solve  (4.1)  and 
(4.2)  by  numerically  integrating  these  equations  along  the  z-axis;  the  details  of  the  method 
have  been  reported  in  chapter  2. 

A  two-dimensional  grid  is  used  by  the  FD  method  to  contain  the  beam  profiles.  In  this  way 
it  is  possible  to  study  the  propagation  of  Gaussian-Hermite  (G-H)  and  Gaussian-Laguerre 
(G-L)  beams  which  are  not  radially  symmetric.  The  field  of  the  G-H  TEMmn  mode,  without 
phase  component,  is  given  by  [34] 


£mn[x>  V>  *)  —  Hm 


(4.3) 
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where  £<>  “  an  amplitude  factor;  Hn{x)  is  the  Hermite  polynomial  of  the  order  n;  and  tu(z) 
is  the  beam  radius  of  the  corresponding  fundamental  Gaussian  beam.  The  G-L  TEMP* 
modes  are  of  the  form  [34] 


**•*■>- sfo 


(4.4) 


where  L*(z)  is  the  Laguerre  polynomial. 

The  M3  factor  of  the  G-H  and  the  G-L  modes  can  be  calculated  from  the  variances  a3  and 
a3  given  by  [32] 


<rl  =  J J(x- i)3I{x,y)dxdy  and  (4.5) 

°l-  jf  ff(y~  v)5/(*>  V) dx  dV »  (4-6) 


where  I(x,  y)  is  the  intensity  of  the  beam  profile;  f  and  y  are  averages  defined  as 

*  =  jjfj  xl[x,y)dxdy  and  y  -  ~  J  J  yl{z,y)dxdy\  (4.7) 

and  N  is  the  normalization  factor  given  by 

N  =  J  J  I{x,y)dxdy.  (4.8) 

These  integrals  can  be  solved  analytically.  For  the  G-H  modes  the  variances  are 


a\  —  w2(l  +  2m)/4  and  a3  =  tu*(l  +  2n)/4 .  (4.9) 

The  variances  of  the  G-L  modes  for  t  £  1  are 

o*  =  o*  =  w,(l  +  2p  +  <)/4;  (4.10) 

for  l  =  1  the  variances  are 


a3  =  3u>*(l  +  p)/4  and  a3  =  w*(l  +  p)/4 . 


(4.11) 
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By  comparing  the  rate  of  change  of  the  variances  with  the  rate  for  a  TEMoo  beam,  the  M * 
factor  is  obtained.  This  gives 

Ml  =  and  M *  =  .  (4.12) 

W2  9  W2 

We  shall  combine  and  to  obtain  M1  using 

A/’  =  MtMt .  (4.13) 


With  this  definition,  M1  is  defined  for  elliptical  beams  and  M1  is  not  changed  by  astigmatic 
lenses. 


The  Gaussian-Schell-model  (GSM)  beam  has  a  Gaussian  intensity  profile  given  by  [35] 


/(r’2)=^)CXP 


(•»(*))  ’ 


(4.14) 


which  is  identical  to  a  fundamental  Gaussian  beam.  However,  the  width  in  (4.14)  differs 
from  the  fundamental  beam  and  is  of  the  form  [36] 


u;2(z)  =  w^l  +  ^  j  with  zR  =  , 


(4.15) 


where  wo  is  the  beam  radius  at  the  waist,  zo  is  the  position  of  the  waist,  and  zr  is  the 
Rayleigh  range.  Equation  (4.15)  shows  that  the  Rayleigh  range  for  the  GSM  beam  is 
reduced  by  the  beam  quality  factor  M2. 


4.2.2  Time  dependence 

In  order  to  compare  the  experimental  results  with  the  numerical  model,  the  time  dependence 
of  the  pump  pulse  has  to  be  taken  into  account.  The  reason  for  this  is  that  the  experimental 
set-up  measures  the  pulse  energy  of  the  Stokes  and  the  pump  pulse  and  not  the  power.  To 
reduce  the  complexity,  the  power  of  the  Stokes  and  the  pump  beam  are  used  instead  of  the 
full  beam  profile.  In  this  way  the  following  equations  are  obtained  [l] 

-dtF  +  d,F  =  A{z)FP, 
c 


(4.16) 
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-dtB  -  dtB  =  A(z)BP ,  (4.17) 

c 

-dtp  +  d.p  =  -^A(z)(F  +  B)P ,  (4.18) 

C  K, 

where  P(z,t),  F(z,t),  and  B(z,t)  are  the  pump  power,  the  forward  scattered  Stokes  power, 
and  the  backward  scattered  Stokes  power,  respectively;  e/n  is  the  speed  of  light  in  the 
medium.  The  gain  parameter  A(z)  is  obtained  from  the  beam  propagation  method,  de¬ 
scribed  above,  using  the  peak  power  of  the  pump  pulse.  To  solve  eqns.  (4.16)-(4.18),  the 
method  of  characteristics  is  used.  The  numerical  method  is  described  in  chapter  2. 

An  experimental  complication  is  the  small  but  significant  reflection  at  the  windows  of  the 
j  Raman  cell.  Because  of  the  very  high  gain  (e*s  at  threshold),  a  small  reflection  can  still 

j  decrease  the  threshold  energy  [22].  Therefore,  the  input  parameters  include  the  reflection 

coefficient  at  the  windows. 

4.3  Results  and  Discussion 

l 

■  4.3.1  Numerical  results  of  the  beam  propagation  method 

We  use  the  beam  propagation  method  of  section  4.2  to  determine  the  amplification  of  a 
Stokes  beam  by  a  nondiffraction-limited  pump  beam.  For  the  initial  Stokes  beam,  a  TEMoo 
beam  is  used.  For  the  pump  beam,  different  shapes  are  used,  namely,  Gaussian- Her  mite 
(G-H)  modes,  Gaussian-Laguerre  (G-L)  modes,  and  Gaussian-Schell-model  (GSM)  beams. 
The  various  pump  profiles  are  given  by  eqns.  (4.3)  and  (4.4),  where  ti/(z)  is  the  beam  radius 
of  the  free  space  Stokes  beam;  the  waist  location  of  the  pump  and  Stokes  beam  is  at  the 
middle  of  the  Raman  cell. 

Figure  4.1  shows  the  amplification  of  the  Stokes  beam  versus  pump  power  for  various 
pump  shapes  with  M7  factors  of  1,  1.73,  3,  and  5.  It  is  clear  that  the  amplifications  near 
the  threshold  level  (e*s)  are  nearly  equal  for  beams  with  different  shapes  but  an  identical 
M7  factor.  For  instance,  the  G-H  TEMoi  mode,  with  an  M1  given  by  (4.13) 

M7  =  MtMv  —  y/l  +  2 m y/TTiii  =>  M7  =  y/l  m  1.73 ,  (4.19) 

has  a  threshold  power  (m  2  a.u.)  almost  identical  to  the  GSM  beam  with  the  same  M7. 


Page 

46 


1 


Figure  4. 1:  Amplification  of  the  Stoke*  beam  vena*  pump  power  for  various  pump  beams.  The 

solid  lines  indicate  the  GSM  beams,  the  dashed  lines  indicate  the  G-H  TEMmn 
beams,  and  the  dotted  lines  indicate  the  G-L  TEMP<  beams. 


By  examining  the  Stokes  beam  profile,  we  see  that  at  the  threshold  level  the  Stokes  beam 
profile  reproduces  to  a  certain  extent  the  shape  of  the  pump  beam.  This  is  illustrated  in 
fig.  4.2  for  a  G-H  TEMoi  pump  beam  of  sufficient  power  so  that  the  amplification  is  at  the 
threshold  level.  It  is  clear  from  the  figure  that  the  Stokes  profile  shows  the  two  maxima  of 
the  pump  profile.  At  low  pump  power,  the  diffraction  dominates  over  the  amplification  so 
that  the  Stokes  beam  does  not  reproduce  the  pump  beam.  This  is  illustrated  in  fig.  4.3, 
where  there  is  only  a  small  stretching  due  to  the  amplification. 

The  relation  between  the  threshold  power  and  the  M1  factor  is,  unfortunately,  not  always 
so  close  as  for  the  pump  shapes  of  fig.  4.1.  In  particular,  the  circular  symmetric  G-L  TEMjo 
mode  shows  a  large  difference  in  amplification  compared  to  the  GSM  beam.  This  is  shown 
in  fig.  4.4  for  the  G-L  TEMio  pump  beam  (p  =  1  and  /  =  0)  with  M*  =  3.  Figure  4.4 
shows  a  much  higher  amplification  for  the  TEMio  beam  compared  to  the  GSM  beam,  which 
results  in  a  threshold  power  that  is  about  40%  lower  than  that  of  the  GSM  beam. 

The  Stokes  beam  profile  for  a  high  power  G-L  TEMio  pump  beam  is  shown  in  fig.  4.5.  It  is 
clear  that  the  Stokes  beam  profile  does  not  reproduce  that  of  the  pump  beam.  The  center 
of  the  pump  profile  is  much  higher  than  the  surrounding  ring  so  that  the  amplification  is 
mainly  at  the  center.  This  particular  profile  combines  a  high  beam  radius  due  to  the  ring 


Figure  4.2:  Beam  profile*,  for  high  pomp  power,  of  the  G-H  TEMoi  pump  beam  and  of  the 

Stoke*  beam  »t  the  entrance,  middle,  and  exit  of  the  Raman  cell. 
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Beam  profiles,  for  low  pump  power,  of  the  G-H  TEMoi  pump  beam  and  of  the 
Stokes  beam  at  the  entrance,  middle,  and  exit  of  the  Raman  cell. 


Figure  4.3: 
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Figure  4.4:  Amplification  of  the  Stokes  beam  versus  pump  power  for  a  G-L  circular  symmetric 

pump  beam  and  a  GSM  pump  beam,  both  with  an  M2  of  3.  The  solid  line  indicates 
the  GSM  beam  and  the  dotted  line  indicates  the  G-L  TEMio  beam. 

with  a  high  intensity  due  to  the  narrow  central  peak;  in  contrast  with  the  Gaussian  shape 
where  a  comparable  high  beam  radius  implicates  a  lower  intensity.  Furthermore,  the  shape 
of  the  profile  does  not  change  while  propagating  so  the  high  amplification  of  the  center  is 
present  during  the  entire  propagation  through  the  Raman  cell.  In  general,  the  shape  of 
an  arbitrary  beam  changes  while  propagating  thereby  reducing  the  effect  of  high  intensity 
peaks. 

An  annular  beam  profile  combines  a  high  beam  radius  with  a  low  intensity  thus  a  low 
amplification  can  be  expected  meaning  a  high  threshold  power.  Such  a  beam  profile  is  the 
“donut”  shape  given  by 

t  (*,  y,  *)  =  £oi(*.  y,  *)  +  it  io(*.  y,  *) ,  (4-20) 

where  £01  and  Eiq  are  the  G-H  modes  of  eqn.  (4.3).  Figure  4.6  shows  the  beam  profiles 
of  the  pump  and  the  Stokes  beam  at  threshold  power.  We  found  a  threshold  power  that 
was  40%  higher  compared  to  a  GSM  beam  with  the  same  Af1  of  2.0.  Since  an  arbitrary 
annular  beam  is  in  general  not  shape  invariant  while  propagating  through  the  Raman  cell, 
the  amplification  of  an  arbitrary  annular  beam  will  be  higher  because  the  low  intensity  ring 
structure  will  be  lost. 
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Figure  4.5:  Beam  profiles,  for  high  pump  power,  of  the  circular  symmetric  G-L  TEMio  pump 

beam  and  of  the  Stokes  beam  at  the  entrance,  middle,  and  exit  of  the  Raman  cell. 
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Figure  4.C:  Beam  profiles  at  threshold  power  of  the  “donut*  pump  beam  and  of  the  Stokes 

beam  at  the  entrance,  middle,  and  exit  of  the  Raman  cell. 
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Figure  4.7:  Threshold  power  versus  M 3  factor.  The  closed  circles  indicate  pump  beams  that 

are  a  random  superposition  of  G-H  modes;  the  solid  line  indicates  the  GSM  pump 
beams;  and  the  open  circles  indicate  the  "donut*  shape  and  the  G-L  TEMio  mode. 

A  realistic  nondiffraction-limited  laser  beam  does  not  consist  of  a  single  high  order  mode 
but  of  a  superposition  of  several  high  order  modes.  Using  a  random  generator,  we  created 
several  superpositions  of  G-H  modes  and  used  the  beam  propagation  method  to  obtain  the 
M2  and  the  threshold  power.  To  construct  the  superposition,  we  used  eqn.(  4.3)  and 

£(*>y»*)  =  y!  cmn£mn(*>  V>*)>  (4.2 1) 

0<m<2  0<n<2 


where  the  complex  phase  factors  Cn»i  are  set  randomly.  Figure  4.7  shows  the  threshold 
power  versus  M2  for  these  superpositions  and  for  the  GSM  beam.  It  is  easily  seen  that  the 
superpositions  have  a  lower  threshold  than  the  GSM  beam.  This  lower  threshold  is  due  to 
the  fact  that  the  superpositions  have  sharp  intensity  peaks  in  the  beam  profile  in  contrast 
with  the  smooth  distribution  of  the  GSM  beam.  Sharp  intensity  peaks  reduce  the  threshold 
power  due  to  their  higher  intensity  and  due  to  the  gain  focusing  effect  in  sharp  peaks  [33,6], 
However,  this  threshold  reduction  is  at  the  most  40%  of  the  GSM  threshold,  therefore  the 
M*  factor  is  still  a  good  measure  for  predicting  the  experimental  threshold.  Additionally, 
figure  4.7  indicates  the  "donut”  shape  and  the  G-L  TEMio  mode  which  shows  that  these 
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Figure  4.8:  Pump  beam  profiles  at  the  entrance  of  the  Raman  cell  with  corresponding  M2 

and  threshold  power  (in  arbitrary  units). 

shapes  are  extreme  cases  due  to  their  particular  profiles. 

The  relation  between  M2  and  the  threshold  power  also  holds  for  profiles  with  a  complicated 
structure.  Figure  4.8  shows  four  pump  profiles  of  a  complicated  structure  which  are  a 
combination  of  G-H  modes  with  a  maximum  order  of  four.  In  the  figure  the  corresponding 
M2  and  threshold  power  (in  arbitrary  units)  are  also  depicted.  The  threshold  powers  are 
lower  than  the  corresponding  thresholds  for  a  GSM  beam  which  are  6.0  and  6.9  (in  a.u.) 
for  the  M2  values  of  6.0  and  7.0,  respectively.  As  explained  in  the  previous  paragraph,  this 
threshold  reduction  is  caused  by  the  sharp  intensity  peaks.  Again  this  reduction  is  within 
40%  of  the  GSM  threshold. 

4.3.2  Experimental  threshold 

Using  the  experimental  set-up  described  in  chapter  3,  we  measured  the  Raman  threshold 
for  two  M2  values  of  the  pump  beam.  The  M2  value  of  the  pump  beam  was  varied  by 
changing  the  diaphragm  in  the  laser.  For  a  diaphragm  of  1.5  mm  diameter,  an  M2  of  1.12 
was  obtained.  An  M2  of  8.7  was  obtained  for  the  6.5  mm  diaphragm.  Figure  4.9  shows  the 


Figure  4.9:  Stokes  pulse  energy  versus  pump  pulse  energy  near  threshold.  The  circles  indicate 

the  experimental  data;  the  solid  lines  indicate  the  computer  simulation. 

threshold  measurement  for  the  two  diaphragms;  as  expected,  the  threshold  energy  increases 
with  M*. 

For  a  correct  computer  simulation,  the  reflection  coefficient  at  the  windows  has  to  be  es¬ 
timated.  Because  of  the  very  high  gain  (»  e18),  small  reflections  of  the  Stokes  beam  at 
the  cell  windows  lower  the  threshold  power  since  the  second  amplification  more  than  com¬ 
pensates  the  reflection  loss  [22].  The  reflection  of  the  baclrward-scattered-Stokes  wavefront 
at  the  entrance  window  is  shown  in  fig.  4.10.  A  fraction  of  the  beam,  typically  0.5%,  is 
reflected  at  the  window,  and  only  a  small  part  of  the  reflected  wavefront  is  amplified  at  the 
focal  region.  The  reflection  coefficient  is  estimated  by  multiplying  the  reflection  coefficient 
of  the  window  (0.5%)  times  the  fraction  of  the  wavefront  that  is  amplified.  This  wavefront 
fraction  is  determined  with  the  beam  propagation  method  from  the  amplification  ratio  of  a 
divergent  and  convergent  Stokes  beam.  In  this  way  the  following  reflection  coefficients  are 
determined:  10-8  for  the  pump  beam  with  Af*  =  1.12,  and  3  x  10-6  for  M5  =  8.7. 

With  the  numerical  model,  we  have  simulated  the  threshold  measurements  of  fig.  4.9.  The 
following  input  parameters  are  used:  the  Raman  gun  coefficient  as  obtained  from  ref.  [26], 
see  chapter  3;  the  reflection  coefficient  at  the  cell  window;  the  measured  pump  power  versus 
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Figure  4.10:  Reflection  of  the  Stokes  wavefront  at  the  entrance  window  of  the  Raman  cell. 

time;  and  a  GSM  pump  beam  with  measured  M*  factor.  Figure  4.9  shows  the  results  of  the 
computer  simulation  as  solid  lines.  It  is  clear  that  the  threshold  energy  of  the  simulations 
match  the  experimental  results.  The  efficiency  above  threshold,  however,  is  not  correct. 
This  is  due  to  other  effects  not  included  in  the  model  that  occur  at  high  pump  energy.  In 
chapter  5  these  effects,  which  reduce  the  conversion  efficiency,  are  studied. 

4.4  Conclusions 

We  have  shown  that  the  threshold  power  of  stimulated  Raman  scattering,  for  a  tightly  fo¬ 
cused  beam,  depends  mostly  on  the  Af*  factor  and  to  a  lesser  extent  on  the  exact  beam  pro¬ 
file.  This  result  has  been  obtained  for  high  order  G aussi&n- Hermite  and  Gaussian-Laguene 
beams;  for  Gaussian-Schell-model  beams;  and  for  several  superpositions  of  Gaussian- Hermite 
modes.  For  a  circular  symmetric  G-L  mode  and  for  a  ring  shaped  G-H  superposition,  the 
largest  difference  in  threshold  power  with  the  GSM  beams  occurs.  Since  this  difference  is 
only  about  40%,  the  Af*  factor  is  a  useful  parameter  for  analysing  the  experimental  thresh¬ 
old.  Furthermore,  for  an  arbitrary  beam  consisting  of  a  superposition  of  high  order  modes, 
we  have  shown  that  the  difference  in  threshold  power  for  beams  with  the  same  Af1  is  often 
smaller  than  40%. 
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The  M2  factor  of  a  Nd:YAG  laser  beam  has  been  measured  in  order  to  analyze  the  SRS 
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threshold  for  a  Raman  cell  containing  high  pressure  methane.  These  experimental  thresh¬ 
olds  correspond  with  the  simulations  using  a  GSM  pump  beam  with  the  measured  M1  values 
of  1.12  and  8.7,  which  confirms  the  usefulness  of  M*  in  analysing  threshold  measurements. 
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5  OPTIMIZATION  OF  CONVERSION  EFFICIENCY 

5.1  Introduction 

Wavelength  conversion  at  high  pump  energies  is  optimized  by  the  use  of  an  astigmatic  lens 
or  a  pair  of  cylindrical  lenses  to  focus  the  pump  beam  into  the  Raman  cell.  Focused  beams 
are  used  to  obtain  a  high  light  intensity,  which  is  required  for  efficient  conversion.  Approxi¬ 
mately,  the  conversion  does  not  change  by  varying  the  focal  length  of  an  ordinary  lens  since 
the  length  of  the  focus  times  the  intensity  at  the  focus  is  constant  [9].  Thus,  a  shorter 
focal  length  results  in  a  shorter  Rayleigh  range  (length  of  focus)  which  is  compensated  by 
the  increased  intensity  at  the  focus.  However,  an  astigmatic  focus  can  be  used  to  modify 
the  gain  and  therefore  the  conversion  [33].  An  astigmatic  focus  gives  a  tuning  facility  by 
varying  A F,  i.e.  the  distance  between  the  two  foci. 

In  the  experimental  study,  a  Nd:  YAG  laser  was  focused  into  a  Raman  cell  containing  88  bar 
of  methane.  The  conversion  of  the  pump  beam  to  other  wavelengths  was  studied  as  a 
function  of  pump  energy  for  an  ordinary  and  for  an  astigmatic  focus.  We  monitored  the 
energy  at  the  following  generated  wavelengths:  the  first  Stokes  at  1.5  fim,  the  second  Stokes 
at  2.8  n m,  the  anti-Stokes  at  0.8  /xm,  and  the  Brillouin  scattering  at  1.06  /xm.  A  dichroic 
mirror,  transparent  for  the  pump  and  highly  reflective  for  the  first  Stokes,  was  placed  at 
the  entrance  of  the  cell  to  enhance  the  first-Stokes  generation. 

Numerical  simulations  were  performed  and  compared  with  the  experimental  results.  All 
the  input  parameters  were  obtained  from  the  literature  and  none  of  the  parameters  were 
fitted  to  match  the  simulations  to  the  experiments.  Ground-state  depletion,  i.e.  reduction 
of  the  density  of  atoms  or  molecules  in  the  ground  state,  was  investigated  numerically.  This 
saturation  effect  occurs  when  the  relaxation  time  of  the  final  Raman  level  is  long  compared 
to  the  pulse  width  of  the  laser,  as  is  the  case  in  our  experiments.  Under  these  conditions, 
the  number  of  generated  Stokes  photons  is  limited  by  the  number  of  molecules  in  the  volume 
where  the  conversion  takes  place. 
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5.2  Theory 


5.2.1  Raman  amplification 


To  reduce  the  complexity,  the  powers  of  the  Stokes  and  the  pump  beam  are  used  instead  of 
the  full  beam  profile.  The  transverse  intensity  is  accounted  for  in  the  model  through  the  area 
of  the  pump  beam  as  a  function  of  z.  The  Raman  propagation  is  then  described  by  coupled 
nonlinear  equations  for  the  pump  power  and  for  the  forward-  and  backward-traveling  Stokes 
powers  [1,37]: 

-dtF  +  d,F  =  A{z){F  +  6)P,  (5.1) 

c 

-d,B- dtB  =  A(z){B  +  6)P,  (5.2) 

c 

-dtP  +  d,P=  -^A{z)(F+  B  +  26)P,  (5.3) 

c  fc, 

where  P(z,  t),  F(z,  t),  and  B(z,  t )  are  the  pump  power,  the  forward  scattered  Stokes  power, 
and  the  backward  scattered  Stokes  power,  respectively;  c/n  is  the  speed  of  light  in  the 
medium;  and  8  is  the  spontaneous  Raman  noise  power.  The  gain  A(z)  is  given  by 

=  <5"> 

where  g  is  the  Raman  gain  coefficient,  A(z)  is  the  cross-sectional  area  of  the  pump  beam  as 
a  function  of  z,  and  ij  is  a  factor  describing  the  gain  focusing  effect.  A  value  of  1.6  for  q  was 
obtained  with  the  beam  propagation  method  described  in  chapter  2.  The  cross-sectional 
area  of  the  pump  beam  is 


A{z)  =  *ws(z)wv(z) , 


(5.5) 


where  2wz  and  2u>v  are  the  beam  diameters  in  the  z-  and  y-direction,  referring  to  the  1/e2 
intensity  points. 

Second-Stokes  generation  by  a  cascade  process,  where  the  first  Stokes  is  shifted  by  SRS, 
can  be  included  by  using  similar  equations  with  the  forward  first-Stokes  power  instead  of 
the  pump  power: 


c 


(5.6) 


TNO  report 


Page 

57 


jdtB,  -  d.Bt  =  Aa(z)(Bt  +  S)F ,  (5.7) 

where  Fa(z,  t)  and  Ba{z,  t)  are  the  forward  scattered  second-Stokes  power  and  the  backward 
scattered  second-Stokes  power,  respectively.  The  gain  Aa(z)  is  given  by  (5.4)  using  the 
Raman  gain  coefficient  g  at  the  first-Stokes  wavelength.  The  propagation  of  the  forward 
first-Stokes  beam  has  to  be  modified  with  depletion  due  to  second-Stokes  generation: 

” d,F  +  d,F  =  A(*)(F  +  $)P 

-±-Aa(z)(Fa  +  Ba  +  26)F.  (5.8) 

Backward  first  Stokes  can  be  ignored  as  a  pump  source  since  the  dichroic  mirror  at  the 
entrance  window  causes  the  forward  first  Stokes  to  be  much  more  intense.  This  assumption 
was  also  checked  numerically;  adding  a  similar  set  of  equations  with  B  as  pump  had  a 
negligible  influence  on  second-Stokes  generation. 

Second-Stokes  generation  by  four-wave  mixing  is  not  taken  into  account  since  the  wave- 
vector  mismatch  is  large  for  high  pressure  methane.  Using  the  data  presented  in  ref.  [38], 
we  obtain  a  wave-vector  mismatch  AJ fc  of  14.6  cm-1  for  CH4  at  88  bar.  The  phase-match 
angle  given  by  yjkj,£k.k/k,ak$  is  30  mrad.  The  divergence  of  the  pump  beam  is  25  mrad 
given  the  beam  radius  of  0.25  cm  at  the  entrance  window  and  a  focal  length  of  10  cm. 
Since  the  divergence  of  the  pump  beam  is  smaller  than  the  phase-match  angle,  four-wave 
mixing  is  suppressed.  Furthermore,  the  pump  beam  is  far  from  diffraction  limited,  which 
also  reduces  four-wave  mixing. 

An  experimental  complication  are  the  small  but  significant  reflections  at  the  windows  of 
the  Raman  cell.  Because  of  the  very  high  gain  (e3s  at  threshold),  a  small  reflection  can 
still  decrease  the  threshold  energy  [22].  Here  the  term  threshold  means  the  pump  power 
corresponding  to  1  percent  conversion  efficiency  [27].  Therefore,  the  input  parameters 
include  the  reflection  coefficient  of  the  windows. 


5.2.2  Ground-state  depletion 

Stimulated  Raman  scattering  induces  a  transition  of  the  methane  molecule  from  the  ground 
state  to  an  excited  state.  If  the  lifetime  of  the  vibrationally  excited  state  is  much  longer 
than  the  pulse  duration,  the  number  of  molecules  in  the  ground  state  will  diminish  [19]. 
Since  each  excited  molecule  corresponds  to  one  scattered  pump  photon  producing  a  single 


Stokes  photon,  the  following  conservation  law  applies: 

dtN+±-dt(IF-IB)  =  0,  (5.9) 

where  N  is  the  number  density  of  methane  molecules  in  the  ground  state,  Ip  is  the  forward 
Stokes  intensity,  Ig  is  the  backward  Stokes  intensity,  and  hu,  is  the  energy  of  the  Stokes 
photon. 

Ground-state  depletion  will  reduce  the  amplification,  which  is  proportional  to  N  and  ex¬ 
pressed  by 


A{z,t) 


JTj_  N(z,t) 
A(z)  '  No 


(5.10) 


where  Nq  is  the  initial  undepleted  density.  Combining  (5.9)  and  (5.10),  we  obtain  for  A: 


dtA  = 


9V 

N0A2hw. 


(dtB  -  dtF) 


(5.11) 


This  equation  is  implemented  in  the  numerical  method  to  study  ground-state  depletion. 


5.3  Experiment 

Details  of  the  experimental  setup  are  described  in  chapter  3.  A  Q-switched  Nd:YAG  laser 
was  used  to  generate  the  pump  beam.  The  laser  provided  1.06  pm  pulses  up  to  350  mJ 
in  energy  with  a  FWHM  pulse  width  of  7  ns.  The  beam  quality,  i.e.  the  M2  factor,  was 
measured  using  a  standard  knife-edge  technique  giving  an  M2  of  8.7.  Details  of  the  M2 
measurement  are  described  in  chapter  3.  The  energy  of  the  laser  beam  was  adjusted  by  an 
attenuator  before  the  beam  entered  the  Raman  cell. 

The  20  cm  long  Raman  cell  contained  methane  with  a  pressure  of  88  bar.  At  the  entrance 
of  the  Raman  cell  a  dichroic  mirror  was  placed  to  reflect  the  backward  scattered  Stokes  into 
the  Raman  cell  in  order  to  obtain  all  Stokes  radiation  in  the  forward  direction. 

An  estimation  of  the  spontaneous  Raman  noise  is  made  using  the  spontaneous  intensity 
rhu,/4A  as  derived  in  [30].  A  noise  power  6  of  10-9  W  is  obtained  in  this  way.  Correcting 
for  the  higher  divergence  of  a  nondiffraction-limited  pump  beam,  we  obtain  a  noise  power 
6  of  M 4  x  10-9  W,  i.e.  10-7  W.  These  literature  values  were  used  in  the  numerical  model 
described  above. 
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Figure  5.1:  Intensity  of  the  pump  on  the  optic*]  axis  for  an  ordinary  focus  (dashed  line)  and 

for  an  astigmatic  focus  (solid  line). 

The  pump  beam  was  focused  into  the  Raman  cell  by  an  ordinary  spherical  lens  of  10  cm 
focal  length  or  by  a  pair  of  cylindrical  lenses  of  10  cm  focal  length  placed  1.3  cm  apart.  For 
a  nondiffraction-limited  beam,  the  Rayleigh  range  is  given  by  zr  —  Jbu>^/(2MJ),  where  tu0 
is  the  beam  radius  at  the  focus.  The  divergence  of  the  focused  beam  is  0.025  rad  given  the 
beam  radius  of  0.25  cm  at  the  entrance  window  and  a  focal  length  of  10  cm.  We  obtain 
a  Rayleigh  range  of  0.47  cm  using  the  expression  for  the  Rayleigh  range  and  the  following 
expression  for  the  divergence:  <(>  =  wq/zr.  With  this  Rayleigh  range,  we  can  calculate  the 
intensity  on  the  optical  axis  given  by 

w2 

/»»(*)  =  h—, - /0-r°  ; - r-=-7Tr  (5.12) 

w  w(z  +  &F/2)w(z  -  AF/2)  ' 

using  the  following  expression  for  the  beam  radius: 

«*(*)  =  «g(i +**/4)-  (513) 

Figure  5.1  shows  the  intensity  on  the  optical  axis  for  the  ordinary  and  the  astigmatic  focus. 
It  is  clear  that  the  intensity  is  significantly  reduced  for  the  astigmatic  focus. 
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Figure  5.2:  Measured  Stokes  energy  versus  pump  energy  for  an  ordinary  focus  and  an  astig¬ 

matic  focus.  The  lines  are  a  guide  to  the  eye  only. 

5.4  Results  and  Discussion 

Using  the  experimental  setup  described  above,  we  measured  the  Stokes  energy  as  a  function 
of  pump  energy  for  a  spherical  lens  and  for  an  astigmatic  lens.  Figure  5.2  shows  the 
experimental  results.  The  maximum  photon  conversion  efficiency  in  the  cell  is  55%  at 
170  mJ  giving  54  mJ  of  Stokes  energy.  Here  we  have  taken  into  account  reflection  losses  of 
5%  at  the  entrance  and  12%  at  the  exit.  The  high  reflection  losses  are  caused  by  the  dichroic 
mirror  at  the  entrance  and  by  the  CaFj  window  and  collimating  lens  at  the  exit.  Uncoated 
CaFj  was  used  in  order  to  measure  the  second  Stokes  at  2.8  pm  which  is  strongly  absorbed 
by  glass.  At  high  pump  energy  an  astigmatic  focus  gives  a  higher  Stokes  energy  compared 
to  an  ordinary  focus.  The  cross-over  energy,  i.e.  the  pump  energy  where  the  ordinary  and 
the  astigmatic  focus  give  the  same  Stokes  energy,  is  104  mJ  as  obtained  from  figure  5.2. 

Radiation  from  competing,  nonlinear  processes  have  also  been  measured.  Figure  5.S  shows 
the  stimulated  Brillouin  scattering  (SBS)  energy.  It  is  clear  that  the  SBS  is  negligible 
since  the  maximum  energy  is  only  2.5  mJ.  Anti-Stokes  radiation  at  0.8  pm  is  shown  in 
figure  5.4  and  is  clearly  negligible.  Second  Stokes  at  2.8  pm  is  not  negligible  at  high  energy 
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Figure  5.4:  Anti-Stokes  energy  venue  pump  energy  for  an  aetigmatic  focus. 
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Figure  5.5:  S«cond-Stoke»  energy  venue  pump  energy  for  nn  astigmatic  focus.  The  solid  line 

indicates  the  simulation. 

as  figure  5.5  shows.  It  has  to  be  noted  that  we  could  only  measure  the  second  Stokes  in 
the  forward  direction  since  the  entrance  window  absorbed  the  2.8  pm.  At  250  mJ  pump 
energy,  the  second  Stokes  is  12  mJ,  which  implies  that  it  is  one  third  in  photon  number  of 
the  first  Stokes  of  70  mJ.  Therefore,  second-Stokes  generation  appears  to  be  an  important 
limiting  factor  for  first-Stokes  generation  at  high  pump  energy. 

Numerical  calculations  were  performed  to  study  the  limiting  processes  for  1.5  pm  generation. 
Because  of  the  very  high  gain  (>  eu),  small  reflections  at  the  cell  windows  lower  the 
threshold  power  since  the  second  amplification  more  than  compensates  the  reflection  loss 
[22].  The  effective  reflections  are  small  since  the  beam  reflected  at  a  window  is  not  focused, 
so  only  the  fraction  of  the  beam  that  is  reflected  into  the  focus  has  to  be  taken  into  account. 
Using  the  procedure  described  in  [39],  we  obtained  the  following  reflection  coefficients:  1.0 
(dichroic  mirror)  and  5  x  10~s  at  1.54  pm  for  the  entrance  and  exit  window,  respectively; 
and  5  x  10~*  at  2.8  pm  for  both  windows. 

Figure  5.6  shows  the  Stokes  energy  as  a  function  of  pump  energy  for  both  lenses  using  the 
numerical  method  presented  earlier.  The  cross-over  energy  according  to  the  simulation  is 
111  mJ,  which  is  in  good  agreement  with  the  experimental  value  of  104  mJ.  Second-Stokes 
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Figure  5.6:  Simulated  Stokes  energy  versus  pump  energy  for  an  ordinary  focus  and  an  astig¬ 

matic  focus. 

simulation  is  depicted  in  figures  5.5  and  5.7,  showing  at  least  a  qualitative  agreement  with 
the  experiments.  The  effect  of  the  astigmatic  focus  is  a  reduction  of  the  gain  (33],  which 
also  reduces  the  gain  for  second-Stokes  generation  resulting  in  a  higher  conversion  above 
the  cross-over  energy. 

The  difference  between  experiment  and  simulation  can  be  attributed  to  the  simplification 
of  ignoring  the  transverse  intensity  of  the  beams.  A  consequence  of  the  transverse  beam 
profile  is  that  the  pump  energy  at  the  beam  edge  is  not  so  easily  converted  as  the  energy  at 
the  center.  This  edge  effect  is  not  taken  into  account  by  the  simulation,  which  results  in  a 
higher  Stokes  energy.  A  more  extensive  numerical  method  using  the  transverse  intensities 
should  give  a  better  correspondence  of  simulation  and  experiment.  We  hope  to  extend  our 
model  in  the  future  with  the  transverse  intensities. 

Ground-state  depletion  is  another  process  that  can  limit  the  conversion  efficiency.  Since 
the  amplification  takes  {dace  predominantly  in  the  focal  region,  a  natural  approach  is  to 
calculate  the  number  of  molecules  in  this  region,  which  should  correspond  to  the  maximum 
number  of  generated  Stokes  photons.  Using  the  Rayleigh  range  x*  and  the  spot  sise  wjf, 
which  were  derived  above,  we  obtain  a  volume  for  the  focal  region  of  4  x  10~4cms  using 
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Figure  5.7:  Second-Stakes  energy  versos  pnmp  energy  for  an  ordinary  foens.  The  solid  line 

indicates  the  simnlation. 

2 zr  x  xwfi.  This  volume  corresponds  to  101*  molecules  giving  a  maximum  Stokes  energy  of 
130  mJ  using  the  Stokes  photon  energy  of  1.3  x  10“ 19  J.  Since  the  generated  Stokes  energy 
is  of  the  same  order  of  magnitude,  we  expected  a  considerable  influence  of  ground-state 
depletion  (GSD).  However,  numerical  calculations  show  a  negligible  influence  of  GSD  when 
the  simulations  with  and  without  GSD  are  compared.  For  an  ordinary  focus  at  a  pump 
energy  of  250  mJ,  a  Stokes  energy  of  133.2  mJ  was  obtained  without  GSD  and  133.1  mJ  was 
obtained  with  GSD  (second-Stokes  generation  was  not  used  in  this  simulation).  Figure  5.8 
shows  (Nq~  W(z))^(z)  versus  z,  which  is  the  quantity  that  indicates  where  the  molecules  are 
excited;  the  area  under  the  curve  is  the  number  of  excited  molecules,  i.e.  101*  corresponding 
with  the  133.1  mJ  Stokes  energy.  It  is  clear  that  the  majority  of  excited  molecules  is  outside 
the  focal  region,  which  results  in  a  much  lower  GSD  effect. 

Experiments  with  several  methane  pressures  were  performed  to  check  whether  the  GSD  is 
negligible.  At  lower  methane  pressures  the  number  of  molecules  that  are  available  for  SRS 
decreases  which  would  result  in  a  lower  Stokes  energy  at  high  pump  energies.  Figure  5.9 
shows  the  Stokes  energy  versus  the  pump  energy  for  methane  pressures  of  50, 70,  and  90  bar. 
It  is  clear  that  there  is  no  significant  reduction  for  low  pressures  at  high  pump  energies, 


z-axis  (cm) 

Figure  5.8:  Ground-itate  depletion  given  by  [No  —  N(*))A  venue  z  (solid  line);  end  number 

density  of  molecules  in  the  ground  state  (dashed  line). 

therefore  GSD  does  not  occur.  The  lower  Stokes  energy  at  low  pump  power  is  due  to  the 
lower  gain  as  a  consequence  of  the  lower  methane  pressure. 


5.5  Conclusions 

It  has  been  shown  that  the  conversion  efficiency  at  high  pump  energy  of  a  Raman  cell  can  be 
improved  by  an  astigmatic  focus.  Our  experiments  show  that  second-Stokes  generation  by 
a  cascade  process  is  the  principal  limiting  factor  for  high  efficiency  first-Stokes  conversion. 
Brillouin  scattering  and  anti-Stokes  generation  were  negligible.  Measurements  at  several 
methane  pressures  indicate  that  the  first-Stokes  generation  is  not  limited  by  ground-state 
depletion.  An  astigmatic  focus  reduces  the  gain,  which  leads  to  a  lower  efficiency  at  low 
pump  power  but  a  higher  efficiency  at  high  pump  power  due  to  the  lower  gain  of  the 
second-Stokes  generation. 

Numerical  simulations  were  performed  and  compared  with  the  experimental  results.  The 
experiments  and  simulations  correspond  reasonably  well  taken  into  account  that  all  input 
parameters  were  obtained  from  the  literature  and  none  of  the  parameters  were  fitted  to 
match  the  simulation  to  the  experiment.  Good  agreement  was  obtained  for  the  cross-over 
energy,  i.e.  the  pump  energy  where  the  ordinary  and  the  astigmatic  focus  give  the  same 
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Figure  5.9:  Measured  Stokes  energy  versus  pump  energy  for  various  methane  pressures.  The 

lines  are  a  guide  to  the  eye  only. 

Stokes  energy,  namely  104  mJ  for  the  experiment  and  111  mJ  for  the  simulation.  Ground- 
state  depletion  was  also  investigated  numerically.  It  has  been  shown  that  the  majority  of 
the  excited  molecules  is  outside  the  focal  region,  which  causes  the  effect  of  ground-state 
depletion  to  be  negligible. 
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6  CONCLUSIONS 


In  this  report  wavelength  conversion  by  stimulated  Raman  scattering  (SRS)  has  been  inves¬ 
tigated  using  experimental  and  numerical  methods.  The  experimental  setup  consisted  of  a 
Nd-.YAG  laair  and  a  high  pressure  methane  Raman  cell  for  the  conversion  of  the  1.06  pm 
radiation  to  the  eye-safe  1.54  pm  wavelength.  Emphasis  is  on  the  identification  and  model¬ 
ing  of  the  relevant  physical  processes  involved  in  the  wavelength  conversion.  Optimisation 
of  the  setup  towards  practical  application  is  described  in  previous  reports  (3,4].  Both  the 
low  and  the  high  pump  power  regime  have  been  investigated. 

A  three-dimensional  finite  difference  method  has  been  presented  for  the  propagation  of 
beams  focused  by  an  astigmatic  lens.  This  FD  method  is  stable  and  unitary,  which  means 
that  the  power  of  the  beams  is  conserved.  The  FD  method  is  used  for  the  modeling  of  a 
Raman  amplifier.  The  results  of  the  FD  method  have  been  cross-checked  with  analytical 
models  and  with  a  numerical  model  that  describes  gain  focusing.  It  has  been  shown  that 
the  FD  method  gives  reliable  results  both  for  low  and  for  high  pump  power. 

Numerical  simulations  have  shown  that  the  threshold  power  of  stimulated  Raman  scattering, 
for  a  tightly  focused  beam,  depends  mostly  on  the  M2  factor  and  to  a  lesser  extent  on  the 
exact  beam  profile.  This  result  has  been  obtained  for  high  order  Gaussian-Hermite  and 
Gaussian-Laguerre  beams;  for  Gaussian-Scheil-model  beams;  and  for  several  superpositions 
of  Gaussian-Hermite  modes.  For  a  circular  symmetric  G-L  mode  and  for  a  ring  shaped  G-H 
superposition,  the  largest  difference  in  threshold  power  with  the  GSM  beams  occurs.  Since 
this  difference  is  only  about  40%,  the  M2  factor  is  a  useful  parameter  for  analysing  the 
experimental  threshold.  Furthermore,  for  an  arbitrary  beam  consisting  of  a  superposition 
of  high  order  modes,  we  have  shown  that  the  difference  in  threshold  power  for  beams  with 
the  same  M2  is  often  smaller  than  40%.  The  A/2  factor  of  a  Nd:YAG  laser  has  been 
measured  in  order  to  analyse  the  SRS  threshold  for  a  Raman  cell  containing  high  pressure 
methane.  These  experimental  thresholds  correspond  with  the  simulations  using  a  GSM 
pump  beam  with  the  measured  M 2  values  of  1.12  and  8.7,  which  confirms  the  usefulness  of 
M2  in  analyzing  threshold  measurements. 

It  has  been  shown  that  the  conversion  efficiency  at  high  pump  energy  of  a  Raman  cell  can  be 
improved  by  an  ast;gmatic  focus.  Our  experiments  show  that  second-Stokes  generation  by  a 
cascade  process  is  the  principal  limiting  factor  for  high  efficiency  first-Stokes  conversion.  An 
astigmatic  focus  reduces  the  gain,  which  leads  to  a  lower  efficiency  at  low  pump  power  but  a 
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higher  efficiency  at  high  pump  power  due  to  the  lower  gain  of  the  second-Stokes  generation. 
Numerical  simulations  were  performed  and  compared  with  the  experimental  results.  The 
experiments  and  simulations  correspond  reasonably  well  taken  into  account  that  all  input 
parameters  were  obtained  from  the  literature  and  none  of  the  parameters  were  fitted  to 
match  the  simulation  to  the  experiment.  Good  agreement  was  obtained  for  the  cross-over 
energy,  i.e.  the  pump  energy  where  the  ordinary  and  the  astigmatic  focus  give  the  same 
Stokes  energy,  namely  104  mJ  for  the  experiment  and  111  mJ  for  the  simulation.  Ground- 
state  depletion  was  also  investigated  numerically.  It  has  been  shown  that  the  majority  of 
the  excited  molecules  is  outside  the  focal  region,  which  causes  the  effect  of  ground-state 
depletion  to  be  negligible.  Measurements  at  several  methane  pressures  confirm  that  the 
first-Stokes  generation  is  not  limited  by  ground-state  depletion. 

We  believe  that  our  numerical  model  will  be  a  useful  tool  for  analyzing  and  predicting  the 
performance  of  various  conversion  setups  with  variation  in  wavelength,  pulse  width,  pulse 
energy,  M2  value,  and  Raman  medium. 
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RUNNING  THE  PROGRAM 


At  the  time  of  this  writing,  the  program  resides  in  the  directory  /ant/heun4/bin  on  the 
file  system  of  the  CONVEX  C230.  Copy  the  files  raabo  end  small .  Inp  to  your  directory, 
so  you  can  run  the  program  by  typing 


raabo  saall 


Here,  raabo  is  the  program  name  (Raman  Amplification  Method  with  Beam  Options)  and 
saall  stands  for  the  input  file  saall .  Inp.  During  the  execution,  the  program  writes  the 
following  information  to  your  screen: 

raabo  (Raman  Amplification  Method  with  Beaa  Options)  version  1.0;  03.4.10 


z.axis  start  value  is  : 
z.axis  stop  value  is  : 
place  of  x-focus  10 

grid  reduction  at  x-focus 
place  of  y-focus  10 

grid  reduction  at  y-focus 


0 

20 

Rayleigh  length  x-focus 
0.102 

Rayleigh  length  y-focus 
0.102 


1.04 

1.04 


*************  beaa  propagation  aethod  **************** 


pump  power  :  l.S34e6 


z 

step 

widthX 

widthY  dz/kdx'2 

gA-area 

gA/bA  tot _f lux 

0.000 

0.2000 

3. 12e-3 

3. 12e-3 

0.4096 

7.85e-3 

7.958 

30.68 

2.000 

0.2000 

2.52e-3 

2.52e-3 

0.6307 

6. 10e-3 

7.957 

30.68 

4.000 

0.1000 

1.92e-3 

1.92e-3 

0.5435 

2.96e-3 

7.954 

30.68 

0.000 

0.05000 

1.33e-3 

1.33e-3 

0.5624 

1 .43e-3 

7.949 

30.68 

8.000 

0.02500 

7.98e-4 

7.98e-4 

0.7848 

6. 13e-4 

7.944 

30.68 

10.00 

0.01250 

6.07e-4 

6.07e-4 

0.9737 

2.06e-4 

7.966 

30.68 

12.00 

0.02500 

7.98e-4 

7.98e-4 

0.7848 

5. 10e-4 

7.994 

30.68 

14.00 

0.05000 

1 . 33e-3 

1 . 33e-3 

0.5624 

1.42e-3 

7.991 

30.68 

10.00 

0.1000 

1.92e-3 

1.92e-3 

0.5435 

2.96e-3 

7.987 

30.68 

18.00 

0.2000 

2.52e-3 

2.52e-3 

0.6307 

5.08e-3 

7.984 

30.68 

20.00 

0.2000 

3. 12e-3 

3. 12e-3 

0.4096 

7.83e-3 

7.983 

30.68 

********  time  profile  and  plane-wave  approximation  *********** 


pump  pulse  duration  :  1 . B0e-8 
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pump  poise  energy  :  8 .OOe-3 
pump  sextans  power  :  1.634e6 
gain  focusing  factor  :  1.273 

(ainisua)  nuaber  of  iterations  tine  profile  k  plane  wave:  4803 


alterations 

tine 

energy.in.cell 

0 

0.000 

0.000 

800 

3.00e-9 

7.21e-4 

1707 

s.oee-o 

8.58e-4 

2605 

8.99e-9 

2.05e-4 

3503 

1.20e-8 

1.46e-5 

4402 

1 . 50e-8 

3.48e-7 

4893 

1.63e-8 

0.000 

With  this  output,  it  is  possible  to  check  whether  the  simulation  gives  valid  results.  The 
beam  propagation  method  uses  an  adaptable  As  listed  under  step,  which  is  reduced  at  the 
focus  to  keep  Az/Jfc(Az)*  below  a  certain  constant  (here  1.0).  width!  and  wldthY  are  the 
Az  and  Ay  of  the  two-dimensional  grid;  the  area  of  the  grid  is  listed  under  gA*area.  A 
useful  check  is  whether  the  beams  are  well  mapped  on  the  grid,  which  can  be  monitored  by 
the  column  gA/bA  showing  the  ratio  of  the  grid  area  and  the  pump  beam  area.  For  a  good 
mapping,  this  ratio  should  be  around  10.  Another  test  is  the  photon  flux  given  in  arbitrary 
units  by  tot-flux,  which  should  be  constant. 

In  the  time  profile  and  plane-wave  approximation,  general  information  about  the  pump 
pulse  is  given.  The  energy  in  the  cell  is  monitored  and  the  iteration  stops  at  the  point 
where  this  energy  is  negligible.  Units  are  not  listed  in  the  output,  but  this  is  no  real 
draw-back  since  all  dimensions  have  fixed  units:  length  scales  are  in  centimeters,  time  is  in 
seconds,  power  is  in  Watt,  and  energy  is  in  Joule. 

As  was  already  mentioned,  this  output  was  produced  for  an  input  file  snail. inp.  Here 
follows  the  contents  of  snail .  inp: 

Input  file  for  the  Stinulated  Raman  Scattering  progran:  ranbo. 

This  is  an  exanple  input  file  showing  a  minimal  number  of  entries. 

The  values  are  only  listed  as  an  example  and  do  not  reflect  a  realistic 
experinent . 

The  units  between  []  are  consents  only.  Just  as  the  characters  after  /*. 


i 

i 


TNO  report 


Page 

A3 


GENERAL  INPUT  PARAMETERS 

•ud.z.axis  »  20.0  tea]  /*  length  Raman  call  */ 

puap_place_x_focus  ■  10.0  [cm]  /*  position  of  locnn  (in  x-z  plana)  */ 
pump _x_width/2  »  0.06  [cn]  /*  baas  width  (half  dianeter  at  1/a ‘2  points)  •/ 
puap.wavenuaber  ■  5a4  [1/ca]  /*  wavenumber  in  aadina  (not  fraa  spaca)  ♦/ 
raaan_galn_stokes  *  5a -10  [ea/W]  /*  Raman  gain  eoafficiant  */ 


BEAM  PROPAGATION  METHOD  :  uses  2-dia  grid  for  beaa  profile 

beaa.propagation  *  on  /*  on  or  off;  skip  this  method  whan  off  */ 

/*  these  two  entries  are  not  used  if  time .plane .wave  is  on  */ 
pump .power  -  2.50e6  [W]  /*  aore  than  one  value  is  allowed  */ 
stokes.power  ■  1.00e-19  [W]  /*  aore  than  one  value  is  allowed,  but  use  the 
sane  number  as  with  puap.power  */ 

TIME-PROFILE  AND  PLANE-WAVE  APPROXIMATION 

tiae .plane .wave  *  on  /*  on  or  off;  tiae  profile  and  plane  wave  approx.  */ 

pulse.width  -  6e-0  [s]  /*  Used  when  puap_profile.f ile  is  auto;  FWHM  */ 
punp.energy  «  8.0e-3  [J]  /*  scales  the  puap  power  to  this  energy; 
aore  than  one  value  is  allowed  */ 

OUTPUT  SPECIFICATIONS 

out. silent  -  off  /*  off  aeans  output  to  stdout,  on  means  silent  */ 


It  is  clear  that  a  complete  simulation,  with  beam  propagation  and  time  dependence  simu¬ 
lation,  can  be  done  with  a  limited  number  of  input  parameters.  However,  the  total  number 
of  parameters  used  for  this  simulation  is  much  higher,  but  these  other  parameters  are  given 
a  default  value.  In  the  output  file  saall .  otp,  the  default  values  are  listed.  A  list  of  all 
parameters  are  given  in  appendix  B. 

Now  we  have  arrived  at  the  point  where  we  can  examine  the  results  of  the  simulation.  The 
output  of  the  program  is  written  to  file  saall. otp.  Since  this  file  is  rather  long,  parts  of 
the  file  are  replaced  by  two  lines  with  dots: 


[like  this] 


I 
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j  Here  follows  the  contents  of  smell .  otp 

t 

OUTPUT  FILE 

rmabo  (Ruin  Amplification  Method  with  Beam  Options)  version  1.0;  03.4.10 
information  from  the  input  file  : 

(standard  units  :  length  [cm] ,  time  [sec] .  power  [V] .  and  energy  [J]  ) 
end.z.axis  *20.0 


out.silent  ■  off 

****end  of  information  inpnt  file***** 

. . . [SCREEN  OUTPUT] 

. . . [SCREEN  OUTPUT] 

******  Stokes  beam  ****** 


out _ stokes  :  intensity  and  phase  are  on  optical  axis 

Note  that  Stokes  power  is  reset  after  each  step 


z 

power 

intens 

phase 

widthX 

wldthY 

Icorner  Ndx/wX  Ndy/wY 

0.000 

1 . 53e-4 

0.03007 

0.000 

0.06000 

0.06000 

1.0e-12 

2.60 

2.60 

2.000 

1.58e-4 

0.07780 

-0.0480 

0.03701 

0.03701 

1  -6e-13 

2.66 

2.66 

20.00 

1.67e-4 

0.07338 

-6.090 

0.03081 

0.03081 

7.6e-14 

3.14 

3.14 

end  Stokes  beam 


******  pump  bean  ****** 


out .pump 

z 

intensity  and  phase  are  on  optical  axis 
power  intens  phase  widthX  widthY 

Icorner  Ndx/wX  Ndy/wY 

0.000 

1 . 634e6 

3.907e8 

0.000 

0.06000 

0.06000 

0.01864 

2.60 

2.60 

2.000 

1 . 534e6 

6.016e8 

-0.0306 

0.04030 

0.04030 

1.86e-3 

2.60 

2.60 

20.00 

1 . 634e6 

3 . 023e8 

-2.816 

0.04002 

0.04002 

8.92e-4 

2.60 

2.60 

end  pump  beam 


****  Time  Profile  and  Plane  Wave  Approximation  **** 


TNO  report 


Page 

A.5 


tiae  profile:  pup  and  Stoke*  (F  and  B) 
n tuber  of  iterations  per  tine  display:  440.1 


tiae 

P.in 

P_out 

F.out 

B.out 

E.in 

E.out 

0.000 

0.960 

0.000 

0.000 

0.000 

0.000 

0.000 

1.501e-9 

4.437e5 

1.497e5 

7.61e-7 

O.OOe-7 

2.33e-4 

4.20e-6 

2.900e-0 

1 . 238e0 

9.000e5 

0.04132 

0.06403 

1.61e-3 

7.0&e-4 

4.496e-0 

1.630e0 

1 . 492e0 

177.3 

183.6 

3.68e-3 

2.67e-3 

6.904e-0 

1 . 176e0 

1.389e6 

40.34 

32.00 

6.76e-3 

4. 01e-3 

7.492e-9 

0.250e5 

8 . 628e6 

0.02231 

0.01703 

7.10e-3 

6.61e-3 

8.990e-9 

2.414e5 

3 . 823e5 

2.34e-6 

2 . OOe-6 

7.72e-3 

7 . 62e-3 

1.049e-8 

6. 931e4 

1 . 248e& 

4.97e-7 

4.73e-7 

7.93e-3 

7.87e-3 

1 . 199e-8 

1 . 503e4 

3.001e4 

6.62e-8 

6.33e-8 

7.09e-3 

7.07e-3 

1.349e~8 

2476 

6681 

8 . 61e-0 

8 . 25e-0 

8.00e-3 

8.00e-3 

1 . 498e-8 

157.6 

810.1 

1 . 17e-9 

1.14e-9 

8.00e-3 

8.00e-3 

Total  energies  pup  and  Stokes. 

Energy  transaitted  pup  :  7.000e-3 
Energy  forward  Stokes  :  3.784e-7 
Energy  backward  Stokes  :  3.660e-7 

E_P+kP/kS*(E_B+E_F)  :  8.000e-3 

************************************************* 

Used  DEFAULT  values: 

out .debug  ■  off 
out_data_f iles  *  off 


reflection.entrance.stokes  *0.0 
ref lection.exit.stokes  “  0.0 

************************************************* 

Value  of  AUTOMATIC  variables: 

display_step_size_z_axis  ■  2 
x_spacing.gr id.points  *  0.003126 
y_spacing.gr id.points  ■  0.003125 

pup_profile_file  *  no  file;  profile  is:  tiae‘2*exp[-(l . 16*tiae/pulsewidth)*2] 
display _tiae .step  •  1.40813e-00 
spontaneous .scattering  »  0 
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************************************************* 

Information  from  the  Input  file  that  was  NOT  used: 

punp .power  •  2.50e6 
stokes.power  ■  1.00e-19 

The  output  file  shown  above  is  the  shortest  possible,  since  all  output  parameters  are  off 
by  default.  Since  the  screen  output  is  repeated  in  the  output  file,  no  information  is  lost 
by  setting  out_silent«on.  For  this  setting,  screen  output  is  suppressed  and  the  program 
can  run  in  the  background  or  in  a  batch  job.  Parallel  processing  is  possible  by  running  the 
program  for  different  input  files  in  the  background. 

Graphical  output  is  not  provided  by  the  program,  but  setting  out  jdata  Jtiles*on  gives  data 
files  with  extension  dto,  which  can  be  easily  read  by  various  graphical  software  packages 
(e.g.  UNIGRAPH). 
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LIST  OF  INPUT  PARAMETERS 


All  parameters  are  listed  in  the  file  full .  lap.  If  the  program  is  modified,  full .  inp  should 
also  be  modified  to  provide  an  up-to-date  list  of  parameters.  Here  follows  the  contents  of 
full . inp: 


Input  file  for  the  Stimulated  Raaan  Scattering  program:  raabo. 

This  is  an  example  input  file  showing  all  the  possible  entries. 

The  values  are  only  listed  as  an  example  and  do  not  reflect  a  realistic 
experiment . 

The  default  values,  if  available,  are  listed  in  the  comments. 


GENERAL  INPUT  PARAMETERS 

end.z.axis  *  20.0  [cm]  /*  length  Raaan  cell  */ 

pump_place_x_f ocus  ■  10.0  [cm]  /*  position  of  focus  (in  x-z  plane)  */ 
pump_x_width/2  ■  0.0572  [cm]  /*  beam  width  (half  diameter  at  l/e*2  points)  */ 
width_at_f ocus  *  off  /*  on  or  off;  width  is  specified  at  focus  or  entrance; 

default:  off  */ 

pump_place_y_f ocus  *  10.0  [cm]  /*  default:  pump _place_x_f ocus  */ 
pump_y_width/2  *  0.0572  [cm]  /*  default:  pump.x.width/2  */ 

pump .beam .quality  *  1.0  /*  factor  ■  l/M‘2  <■  1  ;  contains  ratio  focus-width  or 
divergence  compared  to  a  diffraction  limited  beam;  default:  1.0  */ 
puap.heruite.x  ■  1  /*  order  Heraite  polynomial;  default:  0  */ 
puap_hermite_y  ■  0  /*  default:  0  */ 

/*  these  two  entries  are  used  if  the  heraite  orders  are  0  and  0  */ 
pump.laguerre .angular  *  0  /*  angular  order  Laguerre  mode;  default:  0  */ 
puap.laguerre.radial  *  0  /*  radial  order  Laguerre  node;  default:  0  */ 

puap.wavenuaber  *  6.165e4  [1/ca]  /*  wavenumber  in  medium  (not  free  space)  */ 
stokes.wavenuaber  “  6.165e4  [1/ca]  /*  wavenumber  in  medium  (not  free  space); 

default:  puap.wavenuaber  */ 


raaan_gain_stokes  ■  4.91e-10  [ca/V] 

randoa.seed  *0.1  /*  seed  for  random  generator:  >0  and  <1;  default:  0.1234  */ 
•average  »  1  /♦  number  of  runs  to  average;  default:  1  */ 
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BEAM  PROPAGATION  METHOD  :  uses  2-dia  grid  lor  beu  profile 

be an .propagation  ■  on  /*  on  or  oil;  skip  this  ssthod  whan  oil  */ 
use.syaaetry  •  on  /*  on  or  oil;  on  incrsssas  spaad;  default:  on  */ 

step.size.z.axis  ■  0.2  [ca]  /*  step  size  ol  FD- method;  default:  anto  •/ 
display.step.size.z.axls  ■  2.0  [ea]  /*  rains  or  anto;  default :  anto  */ 
variable .step  -  oil  /*  on,  oil,  or  anto;  dalanlt:  anto  */ 
var.step .criterion  ■  1.0  /*  need  il  variable.step  is  anto;  Units  step  size; 

delault:  1.0  */ 

/*  these  two  entries  are  nsed  il  variable.step  is  on  */ 
var.step.size  ■  0.1  0.05  0.02  0.02  0.02  0.05  0.1  [ca] 

var_8tep_start_z  “0.0  6.0  7.0  8.0  12.0  13.0  15.0  [ca] 

Ngrid.points.on.axis  >  30  /*  the  sane  lor  the  x-  and  y-axis;  default;  40  */ 

/*  grid  spacing  can  be  obtained  autonatically  il  grid.adjnstaent  is  pump  */ 
x_spacing_grid_points  *  0.0055  [ca]  /*  size  delta  x  or  anto;  delault:  anto  */ 
y.spacing_grid_points  »  0.005  [ca]  /*  size  delta  y  or  anto;  delault;  auto  */ 
grid.adjnstaent  “  puap  /*  none,  puap.  stokes,  or  locus;  delault:  pump  */ 
stokeB_grid  ■  oil  /*  on  or  oil;  separate  curvature  lor  stokes;  delault:  oil  */ 
/*  these  lour  entries  are  used  il  grid.adjnstaent  is  ‘locus*  */ 
grid_place_x_focus  ■  10.0  [ca] 

grid_reduction_x_locu8  »  0.02  /*  smaller  than  one  */ 
grid_place_y .locus  -  10.0  [ca] 

grid.reduction.y_f ocus  ■  0.02  /*  smaller  than  one  •/ 

pump .power  “  2.50e6  [W]  /*  wore  than  one  value  is  allowed  */ 
stokes.power  >  1.00e-19  [W]  /*  aore  than  one  value  is  allowed,  bnt  nse  the 

saae  nnaber  as  with  pnap.power  */ 

stokes.noise  >  oil  /*  on  or  oil;  spontaneous  Stokes  power;  default:  oil  */ 

8tokes_place_x_locus  ■  10.0  [ca]  /*  default:  puap_place_x_locus  */ 
stokes_x_width/2  ■  0.0572  [ca]  /*  delault:  pnap_x_width/2  */ 
stokes.place.y .locus  ■  10.0  [ca]  /*  delault:  puap.place.x_l ocus  */ 
stokes.y.width/2  ■  0.0572  [ca]  /*  delault:  puap_x_width/2  */ 
stokes_heraite_x  ■  0  /*  order  Heraite  polynomial;  delault:  0  */ 
stokes.hermite.y  ■  0  /*  delault:  0  */ 

/*  these  two  entries  are  used  il  the  heraite  orders  are  0  and  0  */ 
stokes .laguerre .angular  ■  0  /*  angular  order  Laquerre  mode;  delault:  0  */ 
stokes. laguerre .radial  ■  1  /*  radial  order  Laquerre  node;  delault:  0  */ 

puap.shape  ■  normal  /*  normal,  donut,  tophat,  saooth.tophat ,  GSM; 

delault:  noraal  */ 

stokes.shape  ■  noraal  /*  noraal,  donut,  tophat,  saooth.tophat,  GSM; 
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default:  noraal  */ 

fixed.power.GSM  *  on  /*  used  if  beam.shape  ia  GSM;  default:  on  *>/ 

•dge .transparency  ■  off  /*  on  or  off;  default:  off  */ 

•dge .absorption  -  0.1  [1/cs]  /*  absorption  coefficient;  default:  0.0  */ 

/*  iedge .points  is  used  when  edge .absorption  is  greater  tban  0.0  */ 
•edge .points  ■  1  /*  noraal ly  one,  greater  or  equal  to  one;  default:  1  */ 
low_paas.fi Iter  »  off  /*  on  or  off;  default:  off  */ 

f ully.explic it  -  0.2  /*  if  dx/kdx‘2  <  fully.expl,  the  f.e.  method  ia  used; 

default:  0.0  */ 


TIME-PROFILE  AND  PLANE-WAVE  APPROXIMATION 

tiae.plane.wave  *  on  /*  on  or  off;  tine  profile  and  plane  wave  approx.  */ 

punp.prof ile.f ile  “  auto  /*  naae  data  file;  default:  auto  */ 

pulse .width  *  6e-9  [a]  /*  Used  when  punp.prof ile.f ile  is  auto;  FWHM  */ 

pump.energy  “  8.0e-3  [J]  /*  scales  the  puap  power  to  this  energy; 

more  than  one  value  is  allowed  */ 

tinterva .s.z.axis  -  200  /*  z-axis  la  divided  in  parts  of  equal  lengths; 

default:  200  */ 

gain.window  -  0.1  /*  0<-factor<l ;  anplif ication  only  for 

gain  >  gain_window*gain_nax;  default:  0.0  */ 
gain .focusing  ■  1.674  /*  factor  for  describing  gain.foc;  default:  1.3  */ 

display .time.step  *  0.26e-9  [sec]  /*  value  or  auto;  default:  auto  */ 

spontaneous.scattering  *  le-13  /*  factor,  works  only  at  entrance; 

Stokes.power  is  lactor  *  pump .power;  default:  auto  */ 
spontaneous .power  »  le-7  [W]  /*  used  when  transient .gain  is  off; 

default:  le-7  [W]  */ 

refractive.index  -  1.0  /*  for  speed  of  light  in  Raman  cell;  default:  1.0  */ 

ref lection .entrance .stokes  «  0.0  /*  fraction  backward  to  forward  Stokes; 

default:  0.0  */ 

ref lection.exit.stokes  »  1.0  /*  fraction  forward  to  backward  Stokes; 

default:  0.0  */ 

transient .gain  ■  on  /*  on  or  off,  with  off  the  steady-state  gain  is  used; 

default:  off  */ 

raman.linewidth  -0.2  [1/cm]  /*  linewidth:  Full  Width  Half  Maximum  */ 


TNO  report 


Page 

B.4 


pump.linewidth  *0.0  [1/ca]  /*  linewidth:  Full  Width  Hall  Maximum  */ 
mode .spacing  *0.06  [1/ca]  /*  spacing  ol  longitudinal  nodes  punp  */ 
langeTin.noise  *  on  /*  on  or  oil;  noise  in  the  Renan  polarization  */ 
correct.lresnel#  *  oil  /*  on  or  oil;  correction  in  noise  lor  Fresnel  nuaber  */ 

ground.state.depletion  -  on  /*  on  or  oil.  with  on  gaingrid  is  used; 

delanlt:  oil  */ 

degeneracy. ratio  •  0.0  /*  ratio,  degeneracy  ground  state  divided  by 

degeneracy  excited  state  */ 

atoa.density  ■  2.53el9  [l/ca*3]  /*  nuaber  density  ol  atoas  or  aolecules  */ 

brillouin.scattering  -  oil  /*  on  or  oil,  with  on  SBS  is  included  */ 
brillouin.gain.stokes  *  23.8e-0  [ca/W]  /*  gain  Stiaulated  Brillonin  Scat.  */ 
brillouin.linewidth  -  0.0014  [1/ca]  /*  linewidth:  Full  Width  Hall  Maxinua  */ 

second.order.ranan  *  oil  /*  on  or  oil ;  on  lor  cascade  second  order  Ranan  */ 
raman_gain_stokes2  »  4.8e-10  [ca/W]  /*  gain  second  order  SRS  */ 
stokes2.wavenumber  *  2.33e4  [1/ca]  /*  wavenumber  in  aediua  (not  Iree  space)  */ 
rellection.entrance_stokes2  *  le-4  /*  Iraction  backward  to  lorward  Stokes  */ 
rellection.exit.8tokes2  *  6e-.  /*  Iraction  lorward  to  backward  Stokes  */ 


OUTPUT  SPECIFICATIONS 

out .silent  *  on  /*  oil  means  output  to  stdout,  on  means  silent; 
delault:  oil  */ 

out.intensity  -  x  /*  oil  or  direction:  x,  y.  or  xy; 

the  intensity  along  direction  is  shown;  delault:  oil  */ 
out .phase  ■  oil  /*  oil  or  direction:  x,  y,  or  xy; 

the  phase  along  direction  is  shown;  delault:  oil  */ 
out_2D_intensity  *  10.0  /*  oil  or  position(s)  on  z-axis;  delault:  oil  */ 
out.wavelront  ■  on  /*  on  or  oil;  M*2,  curvature,  etc;  delault:  oil  */ 

out.amp.grid  *  on  /*  on  or  oil;  shows  amplilication  grid;  delault:  oil  */ 
out_vol_grid  *  on  /*  on  or  o 11;  shows  volume  ol  z-grid;  delault:  oil  */ 
out .power. tiae  ■  oil  /*  oil  or  time  point (s)  to  display  powers  in  cell; 

delault:  oil  */ 

out .polar .time  -  le-9  /*  oil  or  time  point (s)  to  display  polarization; 

delault:  oil  */ 

out.polar.index  *  100  /*  oil  or  index  ol  z-position(s)  to  display 

polarization [z.index]  versus  time;  delault:  oil  */ 


out .data .1 lies  ■  oil  /*  on  or  oil;  writes  inlormation  to  data  liles  with 

extension  .dto;  delault:  oil  */ 
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out_reduce_data  *  auto  /*  oil,  auto,  or  reduction  lactor  (integer); 

reduces  data  in  data  liles;  default:  auto  */ 

out .debug  *  oil  /*  on  or  oil;  shows  debugging  inlornation;  delault:  oil  */ 
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DETAILS  OF  THE  PROGRAM 


The  program  is  written  in  the  language  C  using  the  version  of  Kernighan  and  Ritchie  and 
not  the  new  ANSI-C  version.  Using  K&R-C  means  that  the  program  should  be  portable  to 
all  computer  systems  with  a  C-compiler.  For  instance,  we  have  ported  the  program  to  an 
MS-DOS  PC,  where  the  program  ran  (although  not  so  fast)  with  only  a  few  modifications, 
which  were  mostly  caused  by  the  crippled  memory  management  of  the  PC. 

The  source  code  of  almost  10.000  lines  is  divided  among  14  files  not  counting  the  header 
files.  A  special  header  file  is  ramandef .  h  that  contains  the  constant  definitions: 

/*  Header  file  that  contains  the  definitions  of  the  raman  program.  */ 

•define  L0CAL_VAR  0  /*  0  for  maximum  use  of  global  memory  (saves  memory), 

1  for  maximum  use  of  local  memory  (increases  speed)  */ 

•define  NTRD  121  /*  the  maximum  size  of  the  beam  array  */ 

•define  NZ  801  /*  the  max.  number  z-intervals  in  the  plane  wave  approx.  */ 

•define  PI  3.141592654 

•define  h_Planck  6.6262e-34  /*  J.s  */ 

•define  cJLight  2.9979el0  /*  cm/s  */ 

•define  ONE  1.00001  /*  compensates  round-off  in  comparisons  */ 

•define  SIGNAL  -1.234667  /*  signal  value  */ 

•define  NL0WPASS  2  /*  low_pass_f ilter  after  2  steps  */ 

•define  equaKzl,  z2)  (0NE*(zl)>*(z2)  Mt  0NE*(z2)>*(zl))  /*  equal  */ 


The  14  source  files  are  in  alphabetic  order 


average . c 

contains  functions  for  averaging 

inputf il.c 

functions  for  reading  the  input  file 

seaman . c 

for  memory  management 

outfil.c 

basic  functions  for  writing  to  the  output  file 

raaldia . c 

plane-wave  model,  core  functions 

raaldout.c 

plane-wave  model,  output  functions 

ramanout . c 

output  functions  for  beam  propagation  method 

raacontr .c 

contains  upper  level  functions 

ramf unc . c 

beam  propagation  method,  core  functions 

raafunc2.c 

beam  propagation  method,  core  functions 

ramp . c 

highest  level,  contains  main() 

raadomf . c 

random  generators 

rcoaplex.c 

functions  for  complex  numbers 

wavefron.c 

beam  wavefront  analysis 
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At  the  outermost  level,  the  program  works  according  to  the  following  flow-diagram: 


read  input  file _ 

initialization  of  variables  using  input  file 


while  averaging  or  multiple  pump  powers/energies 


yes 

use  beam  propagation  method?  j 

no 

initialize  pump  and  stokes  beam 

beam  propagation  step 

z  =  z  +  A  z 

while  z  <  length  of  raman  cell 

yes 

use  plane-wave  model?  j 

no 

initialize  fields  (and  polarization) 

yes 

\\  beam  propagation  was  used?  j 

no 

amplification  from  BPM 

calculate  amplification 

plane-wave  time  step 

while  energy  in  cell  >  0 

averaging  was  used? 


yes 

calculate  average  information 


no 
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